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1.0  Preliminaries 

1.1  Definition  of  Scientific  Problem: 

This  report  is  an  investigation  of  time-frequency  representations  of  radar  signals 
with  a  view  to  accurately  detect  small  differences  between  radar  returns.  This  objective  is 
congruent  with  BMDO  objectives  of  being  able  to  detect  target  from  decoy  using  active  or 
passive  sensing  methods. 

1.2  The  Aim 

The  aim  of  this  work  is  to  provide  computational  methods  in  signal  analysis  which 
will  permit: 

•  the  detection  of  both  target  scattering  centers  and  resonances. 

•  ultrahigh  resolution  imaging  of  targets. 

•  the  use  of  pattern  recognition  techniques  in  resolving  small  differences  in  surface 

conduction  between  targets  and  decoys. 

1.3  State  of  the  Art  Benchmark  from  which  BSEI  Research  Proceeds 

Radar  means  ranging  and  detection  and  the  detection  process  is  a  level  detection. 
The  signal  used  to  “paint”  a  target  is  usually  of  a  duration  that  in  length  it  is  longer  than  the 
target  length.  Therefore  the  target  becomes  a  “point  scatterer”  and  level  detection  of  the 
target  is  all  that  can  be  used.  If  the  target  is  stationary  -  which  is  not  the  case  in  the  BMDO 
scenario  -  multiple  signal  returns  can  be  obtained  from  the  target  and  using  cross  range  and 
doppler  analysis  a  synthetic  aperture  radar  (SAR)  can  achieve  target  imaging.  Imaging  can 
also  be  achieved  if  the  radar  is  stationary  and  the  target  is  moving  -  a  process  which  is  the 
geometrical  inverse  of  SAR  and  and  known  as  inverse  synthetic  aperture  radar  or  ISAR. 
BMDO  does  utilize  stationary  radars  but  only  late  in  the  trajectory  of  an  incoming  missile. 
In  both  instances  of  SAR  and  ISAR  the  synthetic  aperture  methods  change  the  target  from 
being  a  point  scatterer  to  a  multiple  scatterer.  Furthermore,  the  conventional  methods  of 
target  imaging,  SAR  and  ISAR,  require  not  merely  a  wide  signal  bandwidth,  but  also  a 
long  imaging  time.  Clearly,  the  BMDO  mission  is  served  optimally,  if  there  are  methods 
for  causing  the  target  to  be  a  multiple  scatterer,  but  which  require  neither  the  radar  nor  the 
target  to  be  stationary,  and  which  also  do  not  require  a  long  imaging  time.  In  fact,  there  is  a 
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such  a  method  and  it  requires  the  transmitted  signal  in  space  to  be  less  than  the  length  of  the 
target. 

If  the  transmitted  signal  in  space  is  less  than  the  length  of  the  target,  the  signal 
return  from  the  target  decomposes  into  its  individual  scattering  components.  That  is  to  say, 
there  is  a  unique  signature  from  the  nose,  the  fins,  the  tail,  etc.  Furthermore,  if  the  signal 
return  is  not  averaged  but  is  received  as  a  fast  sampled  time  series  preserving  amplitude, 
phase  frequency  and  perhaps  polarization  modulations,  i.e.,  is  homodyned,  as  opposed  to 
heterodyned,  the  speed  of  the  currents  set  up  on  the  surface  of  the  target  would  have  a 
unique  signature.  In  other  words,  there  is  reason  to  expect  that  the  surface  material 
composition  would  differentiate  signals  returns  from  targets  identical  in  size  but  differing  in 
material  composition.  In  which  case,  a  signature  would  be  available  distinguishing  a 
missile  from  a  decoy,  which  has  a  different  material  composition.  Also,  the  imaging  time  is 
far  less  than  that  required  by  SAR  and  ISAR  and  could  be  on  the  order  of  one  pulse 
transmit  time. 

Given  a  radar  which  transmits  signals  shorter  than  the  length  of  designated  targets, 
which  -  as  it  performs  more  than  just  ranging  and  detection  -  might  be  called  a  sensor,  the 
next  problem  to  be  addressed  is  that  of  return  signal  analysis.  Because  in  SAR,  ISAR  and 
very  short  pulse  radar  the  target  consists  of  a  number  of  scatterers  -  discontinuities, 
cavities,  comers,  etc.  and  each  provides  a  different  backscattering  behavior,  a  first  set  of 
signature  characteristics  identifying  the  target  is  the  backscattering  (from  a  multiple 
scatterer).  A  second  set  of  signature  characteristics  consists  of  the  resonances  of  the  target. 
To  preserve  in  analysis  the  first  set  of  characateristics  requires  precision  in  the  time 
dimension;  to  preserve  in  analysis  the  second  set  of  characteristics  requires  precision  m  the 
frequency  dimension.  The  state  of  the  art  in  analysis  would  be  to  use  wavelet  methods  for 
precision  in  the  time  dimension  and  Fourier  analysis  for  precision  in  the  frequency 
dimension.  The  present  work’s  objective  is  to  provide  (adaptive)  methods  which  provide 
precision  in  both  the  time  and  frequency  dimensions. 

The  requirement  of  analysis  precision  in  both  time  and  frequency  is  occasioned  by 
other  factors.  A  very  short  pulse  interacting  with  a  larger  target  elicits  there  types  of 
responses:  (1)  the  early  time  response;  (2)  the  resonance  response;  and  (3)  the  late  time 
response.  The  first,  the  early  time  response,  is  due  to  currents  being  set  up  on  the  surface 
of  the  target  and  is  target  aspect-dependent.  The  second,  the  resonance  response,  is 
dependent  on  the  length  of  the  target  and  is  aspect-independent  in  harmonic  components  but 
aspect-dependent  with  respect  to  the  amplitude  of  those  harmonic  components.  The  third, 
the  late  time  response,  is  due  to  the  delayed  “ring  down”  of  the  target  after  “painting”  by  the 
signal  and  its  physics  is  similar  to  that  of  the  resonance  response.  Clearly,  target 
identification  should  be  based  on  a  mix  of  the  three  kinds  of  response  with  their  differemg 
aspect  dependencies. 

Another  major  departure  of  the  present  work  from  the  state  of  the  art  is  in  the  nature 
of  the  method  of  target  identification.  The  method  is  neither  level  detection  as  in 
conventional  radar,  nor  imaging  as  in  conventional  SAR  and  ISAR.  The  new  method  is 
pattern  identification.  Given  time-frequency  methods  preserving  both  the  backscattering 
center  signatures  and  resonances  from  a  multiple  scattering  target,  the  identification  method 
would  be  by  pattern  identification  in  the  time-frequency  plain. 

In  summary,  we  offer  methods  different  from  state-of-the-art  for  solving  die 
BMDO  problem  of  providing  an  ultrahigh  resolution  sensor  which  can  identify  missile 
from  decoy.  We  advocate  pulses  shorter  than  the  length  of  the  target,  adaptive  time- 
frequency  analysis  methods  which  preserve  target  backscattering  centers  and  resonances, 
and  pattern  identification  of  target  and  decoy. 
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1.4  Deficiencies  in  State  of  the  Art  Addressed  by  BSEI  Research 

BMDO  lacks  an  active  high  resolution  sensor  as  well  as  the  capability  of  detecting 
missiles  from  debris  and  decoys.  A  short  pulse  sensor  is  smaller  in  size,  weight  and  power 
than  conventional  radars  and  requires  a  smaller  antenna.  Therefore  a  short  pulse  sensor  can 
be  deployed  either  on  satellites  or  as  a  missile  seeker. 

The  state  of  the  art  in  conventional  radars  can  only  achieve  high  resolution  in  the 
SAR  or  ISAR  configuration.  The  former  requires  a  stationary  target  with  a  moving  radar 
and  the  latter  requires  a  moving  target  with  a  stationary  radar.  A  short  pulse  system  does 
not  have  such  requirements  to  achieve  high  resolution.  Denied  the  SAR  and  ISAR 
configurations,  conventional  radars  achieve  target  identification  by  level  detection. 
Therefore  a  target  missile  cannot  be  distinguished  from  a  decoy.  A  short  pulse  system 
achieves  a  signal  return  which  can  be  analyzed  using  time-frequency  methods  and  target 
can  be  distinguished  from  decoy  by  pattern  matching  techniques. 

1.5  Implementation  by  BMDO 

The  implementation  of  a  short  pulses  sensor/radar  together  with  the  correct  time- 
frequency  approaches  to  processing  is  straightforward  if  a  development  program  is 
initiated.  Alternatively,  simulation  and  modeling  can  take  place  using  the  VDHTB  testbed. 

1.6  What  is  new,  innovative  in  BSEI's  research  results  and  what  are  the 
implications? 

1.6.1  New  methods  for  achieving  discrimination  of  both  frequency- 
dependent  aspects  of  the  target  and  time-dependent  scattering  centers. 

The  Adaptive  Gabor  Representation  (AGR)  was  used  to  analyze  a  radar  signal 
return  in  two  forms:  the  functions  were  either  varying  Q  or  constant  Q.  We  found  the 
constant  Q  form  is  optimum.  Frequency  hopping  signals  as  well  as  monocycle  pulses  were 
analyzed  with  the  AGR  and  the  Combined  Wigner-Ville  Distribution  of  the  AGR 
calculated.  It  was  shown  that  with  all  frequency  components  combined  the  frequency 
components  of  the  signal  can  be  identified,  but  not  die  pulse.  However,  when  the 
Combined  Wigner-Ville  representation  of  the  AGR  is  calculated  using  only  the  high 
frequency  components  of  the  latter,  it  was  seen  that  the  monocycle  can  be  detected.  Thus, 
these  methods  seems  entirely  appropriate  for  detecting  the  frequency-dependent  and  the 
time-dependent  (transient)  scattering  centers  of  targets. 

The  AGR  and  the  Combined  and  Single  Wigner-Ville  Distribution  of  the  AGR  were 
calculated  for  frequency  hopping  signals,  pulses  and  radar  return  signals  The  AGR  requires 
finding  the  optimum  wavelet  center  frequency,  optimum  dwell  and  optimum  modulating 
envelope.  The  hypothesis  under  test  was  that  whereas  the  AGR  has  been  viewed  as  an 

optimization  with  respect  to  three  variables:  ap,  t0  and  f0,  with  ap  setting  the  envelope  size 
independently  of  the  other  two  variables,  use  of  a  constant  Q  Gabor  wavelet  would  permit 

an  efficient  optimization  using  only  one  variable  -  ap  -  with  the  other  signal  parameters 
functions  of  that  one  variable.  It  is  shown  that  the  constant  Q  Gabor  wavelet  with  f0.t0 
=Af.At  =  1/2,  is  optimum  and  that  only  one  variable  is  required  for  optimization. 

1.6.2  New  wavelet  analysis 

A  new  wavelet  analysis  was  developed  based  on  orthogonal  parabolic  cylinder  or 
Weber  functions  of  increasing  level  described  by  the  expansions:  Af.At  =  f0.t0  =  1/2  (2n  + 
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1)  n=0  1,2, . (i.e.,  Constant  Q  Orthogonal  Weber  Function  Wavelets  (WFWs).  Using 

a  frequency  hopping  and  pulse  test  signals,  these  wavelets  are  able  to  detect  both  resonance 
and  pulse  scattering  features.  Using  scattergram  plots  between  three  levels  of  wavelets  and 
two  scale  dilations,  unique  features  of  both  the  test  signal  and  radar  return  signals  can  be 
demonstrated.  This  research  is  unique  in  two  respects:  using  orthogonal  wavelets  of 
different  levels  and  in  the  scatterogram  cross-comparison  of  three  levels  of  analysis,  which 
can  provide  unique  signatures  for  targets.  In  the  case  of  the  test  signal  of  a  frequency 
hopping  signal,  plus  a  pulse,  it  is  seen  that  discernible  cyclic  patterns  occur  which  may  be 
correlated  with  the  number  of  frequencies  in  the  signal  -  four.  In  the  case  of  the  radar  return 
test  signal,  it  is  seen  that  recognizable  helical  patterns  occur,  which  may  provide  a  unique 

signature  for  targets.  ,  .  ,  .  ... 

Using  a  frequency  hopping  and  pulse  test  signals,  WFWs  permit  detection  of  both 
resonance  and  pulse  scattering  features.  The  characteristics  of  these  wavelets  are  based  on 
IIF  filter  principles.  It  is  customary  to  believe  that  a  discrete  time  scaling  (low  pass)  wavelet 
is  first  necessary  to  obtain  the  discrete  time  (high  pass)  wavelet.  Here  we  show  that  the 
expansion  method  permits  an  estimate  of  die  continuous  time  high  pass  wavelet  from  the 
summed  Fourier  spectra  of  the  expansions,  n  =  1,  2,...  °°.  Future  work  will  examine 
whether  repeated  analysis  of  the  original  signal  by  higher  order  wavelets  is  equivalent  to  the 
customary  method  of  using  but  one  scaling  function  and  one  wavelet  with  signal 
decimation. 

1.6.3  Instantaneous  phase  information. 

Instantaneous  phase  information  was  obtained  from  application  of  wavelet  analysis. 

A  detailed  thumbprint  of  a  complex  (target  return)  signal  can  be  obtained  which  can  be  used 
as  data  supporting  a  detailed  wavelet  modulus  description  of  the  same  signal.  A  result  is 
obtained  by  calculating  the  differential  phase  which  also  provides  supporting  data.  There 
are  also  optimum  dilations  for  providing  correlated  measures. 

The  probability  distribution  function  for  the  target  return  signal  processed  with  a 
specific  dilation  of  the  wavelet  and  the  mutual  information  for  two  dilations  also  provides 
supporting  data  to  the  modulus  thumbprint  of  the  signal. 

The  phase  information  and  mutual  information  provide  excellent  supporting  data  to 
that  providing  by  the  wavelet  modulus.  In  combination  a  more  detailed  analysis  of  target 
information  is  possible. 

1.6.4  Target  Differentiation. 

Differentiation  of  two  target  return  signals  was  accomplished  using  (1)  wavelet 
processing  and  (2)  fractal  analysis  of  the  results  of  that  wavelet  analysis.  Log-Log  plots 
(amplitude  versus  filter  scale)  revealed  stable  distributions  (as  opposed  to  asymptotically 
Gaussian  distribution)  behavior  -  increasing  in  the  case  of  one  signal,  and  decreasing  in  the 
case  of  the  other.  This  is  a  fractal  property.  The  computed  fractional  dimension  per  time  bm 
revealed  clear  differences  between  die  two  signals.  The  Fano  factors  for  means  across 
filters  and  for  each  time  bin  were  also  computed.  There  is  a  clear  difference  in  the  Fano 
factor  patterns  generated  permitting  the  differentiation  of  the  signals.  Thus,  fractal  analysis 
appears  to  be  a  promising  approach  to  differentiating  signals  and  will  be  further  explored  in 
future  months. 

1.6.5  Nonlinear  Collective  Exzcitations. 

Analyses  in  terms  of  (a)  higher-order  symmetries;  (b)  differential  forms;  and  (c) 
group  theoretic,  and  formulations  in  electromagnetic  theory  of  the  role  of  the  A-field  in  (a) 
nonlinear  collective  excitations  in  a  2-dimensional  electron  gas;  and  (b)  solitons  in  new 
SAW  devices  and  heterostructures  in  general.  The  method  of  generation  of  modulated 
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TeraHertz  pulses  utilizes  a  plasma  mirror  dependent  on  the  formation  of  soliton 
transmission. 

1.6.6  The  Bootstrap  Method. 

The  Bootstrap  method  was  used  for  estimating  the  accuracy  of  a  paramater 
estimator  taken,  in  the  instances  addressed,  as  radar  signal  returns.  We  demonstrate  -  for 
the  first  time  to  our  knowledge  -  that  the  technique  can  be  applied  to  distinguish  targets.  We 
also  apply  techniques  taken  from  fractal  theory  and  deterministic  chaos  theory  to  show  that 
radar  returns  have  both  a  fractal  nature  and  are  not  time  series  produced  by  a  random 
mechanism,  but  rather  by  a  deterministic  chaotic  mechanism.  We  are  able  to  show  that  the 
fractal  properties  of  the  power  spectra,  the  variance  and  Fano  Factor  distinguish  between 
two  target  returns.  Use  of  embedding  dimensional  analysis  also  reveals  the  deterministic 
chaotic,  as  opposed  to  random,  nature  of  the  radar  returns. 


1.7  Implications  of  this  research 

This  research  provides  time-frequency  methods  for  analyzing  sensor/radar  signal 
returns  to  achieve  ultrahigh  resolution  and  the  possibility  of  distinguishing  target  from 
decoy.  The  research  also  indicates  the  advantage  of  deploying  a  short  pulse  sensor/radar,  in 
which  the  pulse  is  shorter  than  the  length  of  the  target. 

1.8  Future  work 

The  completed  work  addressed  a  number  of  new  research  areas  which  require  much 
more  development  time  than  was  given  in  the  limited  time  available  in  the  present  contract. 
These  areas  include:  Constant  Q  Orthogonal  Weber  Function  Wavelets,  application  of  the 
Bootstap  method,  fractal  analysis  methods  and  the  Adaptive  Gabor  representation.  The 
development  of  these  areas  can  take  place  in  a  simulation  environment. 

1.9  Future  needs 

To  conduct  future  work  in  these  areas  there  are  two  requirements:  (1)  a  source  of 
realistic  signal  returns  from  missiles  and  decoys;  and  (2)  a  testbed  implementation,  e.g.,  the 
VDHTB.  Alternative  to  the  testbed  implementation  using  computer  analysis  of  the  signals, 
the  algorithms  described  in  the  present  work  could  be  implemented  in  ASIC  or  hardware 
designs.  However,  the  preference  would  be  to  firstly  test  the  algorithms  in  a  testbed  and 
then  proceed  to  hardware  implementation. 


2.0  Summary  of  Results  of  Contract  DAAG55-98-C-0044 

These  reported  investigations  of  time-frequency  representations  of  radar  signals 
were  acrried  out  with  a  view  to  accurately  detect  small  differences  between  radar  returns. 
This  objective  is  congruent  with  BMDO  objectives  of  being  able  to  detect  target  from 
decoy. 

2.1  Representative  radar  data  was  analyzed  with  (1)  continuous  Gabor  wavelets;  (2) 
ambiguity  functions;  and  (3)  Wigner-Ville  distributions,  in  order  to  determine  the  optimum 
form  of  analysis  for  the  RF  sensor  returns. 

2.2  We  are  able  to  show  the  superiority  of  the  wavelet  method  and  the  usefulness  of 
Gaussian  differential  wavelets  in  detecting  local  differences  in  the  smoothness  of  the  decay 
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of  transitions  in  the  signal.  The  capability  of  detecting  such  differences  may  be  critical  in 
discriminating  targets  from  decoys.  In  the  case  of  the  spectrogram,  we  showed  that  the 
method  can  pick  out  well  the  onset  and  offset  of  the  various  modulations  in  the  signal. 
However,  due  to  the  fact  that  the  window  is  constant  (even  although  the  modulated 
frequency  of  the  kernel  is  changing),  the  spectrogram  provides  a  biased  representation  of 
the  signal  as  the  sampling  window  excludes  the  influence  of  signal  components  outside  the 
unchanging  windowed  ubox”.  The  spectrogram  is  thus  not  an  optimum  representation  of 
signals  and  future  work  will  demonstrate  how  to  amend  the  spectrogram  method  to  obtain 
wavelet  representations  from  a  similar  commencement. 

2.3  New  methods  were  developed  for  achieving  heightened  discrimination  capability.  The 
action  of  the  Zak  Transform  in  achieving  this  aim  was  investigated  using  a  radar  return 
signal  exemplar  with  variance  of  the  Zak  summing  variable  k.  It  was  shown  that 
discrimination  capability  of  fine  structure  is  lost  as  the  range  of  k  increases.  In  order  to 
investigate  improvements  in  discrimination  capability,  the  Zak  transform  was  performed 
not  on  the  radar  exemplar  but  on  the  radar  exemplar  after  Gabor  analysis  with  kernels  of 
various  sizes.  It  was  shown  that  there  are  distinctive  patterns  for  each  of  the  kernels  used. 
Therefore  a  radar  return  signal  can  be  analyzed  into  a  4-dimensional  pattern  signature  when 
the  results  of  the  analysis  patterns  of  multiple  kernels  are  summed.  This  method  holds 
promise  for  increasing  fine  discrimination  capability. 

2.4  New  methods  were  developed  for  achieving  discrimination  of  both  frequency- 
dependent  aspects  of  the  target  and  time-dependent  scattering  centers.  The  Adaptive  Gabor 
Representation  (AGR)  was  used  to  analyze  a  radar  signal  return  in  two  forms:  the  functions 
were  either  varying  Q  or  constant  Q.  We  found  the  constant  Q  form  appears  optimum. 
Frequency  hopping  signals  as  well  as  monocycle  pulses  were  analyzed  with  the  AGR  and 
the  Combined  Wigner- Ville  Distribution  of  the  AGR  calculated.  It  was  shown  that  with  all 
frequency  components  combined  the  frequency  components  of  the  signal  can  be  identified, 
but  not  the  pulse.  However,  when  the  Combined  Wigner-Ville  representation  of  the  AGR 
is  calculated  using  only  the  high  frequency  components  of  the  latter,  it  was  seen  that  the 
monocycle  can  be  detected.  Thus,  these  methods  seems  entirely  appropriate  for  detecting 
the  frequency-dependent  and  the  time-dependent  (transient)  scattering  centers  of  targets. 

2.5  New  methods  continued  to  be  developed  for  achieving  discrimination  of  both 
frequency-dependent  aspects  of  the  target  and  time-dependent  scattering  centers.  The 
Adaptive  Gabor  Representation  (AGR)  and  the  Combined  and  Single  Wigner-Ville 
Distribution  of  the  AGR  were  calculated  for  frequency  hopping  signals,  pulses  and  radar 
return  signals  The  AGR  requires  finding  the  optimum  wavelet  center  frequency,  optimum 
dwell  and  optimum  modulating  envelope.  The  hypothesis  under  test  was  that  whereas  the 

AGR  has  been  viewed  as  an  optimization  with  respect  to  three  variables:  <xp,  t0  and  f0,  with 

ap  setting  the  envelope  size  independently  of  the  other  two  variables,  use  of  a  constant  Q 

Gabor  wavelet  would  permit  an  efficient  optimization  using  only  one  variable  -  ap  -  with 
the  other  signal  parameters  functions  of  that  one  variable.  It  is  shown  that  the  constant  Q 
Gabor  wavelet  with  f0.t0  =Af.At  =  1/2,  is  optimum  and  that  only  one  variable  is  required 
for  optimization. 

2.6  Investigations  were  initiated  of  wavelets  based  on  orthogonal  parabolic  cylinder  or 
Weber  functions  of  increasing  level  described  by  the  expansions:  Af.At=f0.t0  =  1/2  (In  + 

1),  n  =  0,1,2, .  (i.e..  Constant  Q  Orthogonal  Weber  Function  Wavelets  (CQOWs). 

Using  a  frequency  hopping  and  pulse  test  signals,  these  wavelets  are  able  to  detect  both 
resonance  and  pulse  scattering  features.  Using  scattergram  plots  between  three  levels  of 
wavelets  and  two  scale  dilations,  unique  features  of  both  the  test  signal  and  radar  return 
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signals  can  be  demonstrated.  This  research  is  unique  in  two  respects:  using  orthogonal 
wavelets  of  different  levels  and  in  the  scatterogram  cross-comparison  of  three  levels  of 
analysis,  which  can  provide  unique  signatures  for  targets.  In  the  case  of  the  test  signal  of  a 
frequency  hopping  signal,  plus  a  pulse,  it  is  seen  that  discernible  cyclic  patterns  occur 
which  may  be  correlated  with  the  number  of  frequencies  in  the  signal  -  four.  In  the  case  of 
the  radar  return  test  signal,  it  is  seen  that  recognizable  helical  patterns  occur,  which  may 
provide  a  unique  signature  for  targets. 

2.7  Further  investigation  sof  CQOWs  showed  that  using  a  frequency  hopping  and  pulse 
test  signals,  these  wavelets  are  able  to  detect  both  resonance  and  pulse  scattering  features. 
The  characteristics  of  these  wavelets  are  based  on  IIF  filter  principles.  It  is  customary  to 
believe  that  a  discrete  time  scaling  (low  pass)  wavelet  is  first  necessary  to  obtain  the 
discrete  time  (high  pass)  wavelet.  Here  we  show  that  the  expansion  method  permits  an 
estimate  of  the  continuous  time  high  pass  wavelet  from  the  summed  Fourier  spectra  of  the 
expansions,  n  =  1,  2,...  Future  work  will  examine  whether  repeated  analysis  of  the 
original  signal  by  higher  order  wavelets  is  equivalent  to  the  customary  method  of  using  but 
one  scaling  function  and  one  wavelet  with  signal  decimation. 

2.8  Methods  were  developed  addressing  instantaneous  phase  information  obtainable  from 
application  of  wavelet  analysis.  We  show  that  a  detailed  thumbprint  of  a  complex  (target 
return)  signal  can  be  obtained  which  can  be  used  as  data  supporting  a  detailed  wavelet 
modulus  description  of  the  same  signal.  We  also  show  a  result  from  calculating  the 
differential  phase  which  also  provides  supporting  data.  It  is  also  shown  that  there  are 
optimum  dilations  for  providing  correlated  measures. 

We  also  calculated  the  probability  distribution  function  for  the  target  return  signal 
processed  with  a  specific  dilation  of  the  wavelet  and  the  mutual  information  for  two 
dilations.  This  form  of  analysis  also  provides  supporting  data  to  the  modulus  thumbprint  of 

the  signal.  . ,  „  . 

The  phase  information  and  mutual  information  provide  excellent  supporting  data  to 

that  providing  by  the  wavelet  modulus.  In  combination  a  more  detailed  analysis  of  target 

information  is  possible. 

2.9  Methods  were  developed  addressing  the  differentiation  of  two  target  return  signals 
using  (1)  wavelet  processing  and  (2)  fractal  analysis  of  the  results  of  that  wavelet  analysis. 
Log-Log  plots  (amplitude  versus  filter  scale)  revealed  stable  distributions  (as  opposed  to 
asymptotically  Gaussian  distribution)  behavior  -  increasing  in  the  case  of  one  signal,  and 
decreasing  in  the  case  of  the  other.  This  is  a  fractal  property.  The  computed  fractional 
dimension  per  time  bin  revealed  clear  differences  between  the  two  signals.  The  Fano 
factors  for  means  across  filters  and  for  each  time  bin  were  also  computed.  There  is  a  clear 
difference  in  the  Fano  factor  patterns  generated  permitting  the  differentiation  of  the  signals. 
Thus,  fractal  analysis  appears  to  be  a  promising  approach  to  differentiating  signals  and  will 
be  further  explored  in  fiiture  work. 

2.10  Analytical  methods  were  developed  in  terms  of  (a)  higher-order  symmetries;  (b) 
differential  forms;  and  (c)  group  theoretic,  and  formulations  in  electromagnetic  theory  of 
the  role  of  the  A-field  in  (a)  nonlinear  collective  excitations  in  a  2-dimensional  electron  gas; 
and  (b)  solitons  in  new  SAW  devices  and  heterostructures  in  general.  The  method  of 
generation  of  modulated  TeraHertz  pulses  utilizes  a  plasma  mirror  dependent  on  the 
formation  of  soliton  transmission.  Therefore  an  analysis  was  commenced  of  such  soliton 
transmissions  using  the  various  approaches  described. 

2.11  Analytical  methods  were  developed  based  on  the  Bootstrap  method  for  estimating  the 
accuracy  of  a  paramater  estimator  taken,  in  the  instances  addressed,  as  radar  signal  returns. 
Here,  we  show  -  for  the  first  time  to  our  knowledge  -  that  the  technique  can  be  applied  to 
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distinguish  targets.  We  also  apply  techniques  taken  from  fractal  theory  and  deterministic 
chaos  theory  to  show  that  radar  returns  have  both  a  fractal  nature  and  are  not  time  series 
produced  by  a  random  mechanism,  but  rather  by  a  deterministic  chaotic  mechanism.  We  are 
able  to  show  that  the  fractal  properties  of  the  power  spectra,  the  variance  and  Fano  Factor 
distinguish  between  two  target  returns.  Use  of  embedding  dimensional  analysis  also 
reveals  the  deterministic  chaotic,  as  opposed  to  random,  nature  of  the  radar  returns. 


3.0  Continuous  Gabor  Wavelets 

The  Gabor  wavelet  kernel  is: 

iff  =  Exp[-£/]  x  Cosfy) . 


and  in  general: 


V(t)  =  g(t)Exp[i£2t], 


where  g(t)  is  a  Gaussian  window: 


g(  0  = 


with  4  = 


1 

2c1' 


The  continuous  Gabor  wavelet  transform  is  then  defined  as: 


Two  examples  of  radar  signal  returns  are  shonw  in  Figs.  3.1A  and  3.2A.  The 

Gabor  transforms,  Wf(t,s),  are  shown  in  Figs.  3.1B-D  and  3.2B-D,  respectively.  It  can  be 
seen  that  the  two  signal  returns  can  be  identified  by  3-D  pattern  detection. 


Tgt:l ,Sig:  #0001, AZ:-24. 9899997 fAl:13. *0999994 


Fig.  3.1A.  1st  example  of  radar  signal  return. 
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Tgt:l  ,sig^#0001  ,li:-24 .9899997,11:13.80999994 


Tgt:l  ,$ig:#0001  ,A.z:  -24 . 9899997,11: 13.80999994 


Tgt:l,Sig:#1162,Az:-22.S30000S,JLl:7. 330000114 


Tgt:l,Sig:#1162,AZJ-22.S300006, *1:7. 330000114 


4.  Ambiguity  Function. 

The  ambiguity  function  is: 


Af{  t,cq)  =  J  f^t  +  ^  f^t  -  ^jExp[-iax]dt , 


or,  in  the  frequency  domain: 


A/(t,<o)  =  — J/^v  +  yjf  v-y  Exp[irv]dv, 

and  measures  the  spread  of  the  signal  /  in  time  and  of  /  in  frequency,  i.e.,  the  energy 

concentration.  Figs.  4.1A-C  represents  the  energy  spread  of  the  radar  return  example 
shown  in  Fig.  3.1  A. 


»w.n, x.rrc.-B.nd  wbigait,  functio.  Tst:l.si9:«001,»i!- 


24. 9899997, Al:13. 80 99 9994 


Xdrrow-Bsna  Amp lguity  mncticc  Tgt:l  ,5ig;#oooi  ,Az:-2*.9899997,Al:i  3.80999994 


.  '  ,  ,  ”il* 

'•  "  . . . 


Fig.  4.1  A,  B  &  C.  First  example:  ambiguity  function 


5.0  Wigner-Ville  Distribution. 

The  Wigner  -Ville  distribution  is  a  quadratic  form: 


Pvf ( t,a> )  =  J f^t  +  Exp[-ion]dt 


or 


Pvf(t,co) 


=I;H> 


* 1  (o  -  —  \Exp[-i vt]dv 

Ami  j 


and  defines  a  time-varying  power  spectrum  for  non-stationary  processes.  Simply  stated, 
the  Wigner-Ville  transform  is  the  cross-correlation  of  a  signal  with  itself  after  a  time  and 
frequency  shift. 

The  two  characteristics  which  make  the  Wigner-Ville  distribution  nonoptimal  for 
signal  analysis  purposes  are  that  it  can  take  on  negative  values  and  that  the  presence  of 
interference  terms.  Fig.  5.1  A.  is  an  example  of  two  Gabor  signals.  Figs.  5.1B-D  show  the 
Wigner-Ville  transforms  of  those  signals.  The  interference  terms  are  clearly  seen. 
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5.1B.C&D.  Wi 


Fig.  5.2A-D  are  the  Wigner-Ville  transforms  of  the  radar  return  example  shown  in 
Fig.  3.1  A. 


Vigner-Ville  Function  Tgt:i  ,sig:#oooi  ,*-z:-2«.9899997,Al!l3. 8099999* 
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Higner-vUle  Function  Tgt:l  ,Sig:#oooi  ,Az:-2$. 9S99997,Al:l3 . 8Q999994 


C. 

Fig.  5.2A,  B  &  C.  Wigner-Ville  Transform  of  1st  Example. 
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6.0  Wavelets  compared  with  spectrogram 
6.0.1  Prologue. 

The  radar  signal  shown  in  the  following  Fig.  6.1  is  wavelet  transformed  in  the 
examples  to  follow.  The  frequency  chirp  shown  in  Fig.  6.2  is  used  as  an  example  in  the 
spectrogram  analysis. 


*rgt:l ,Sig:#0001 ,AZ!-2*.9$99997,Al:13. 8099999* 


Fig.  6.1.  Radar  signal  example. 


Fig.  6.2.  Linear  chirp  f(t)  =  Exp[iar2]. 
The  wavelet  integral  is  defined: 


Calderon  and,  independently,  Grossman  and  Morlet  proved  that  if  yr  e  L2(SR)  is  a  real 

function  (that  is,  L2 (SR) represent  finite  energy  functions  such  that:  J|/(r)fdr  <  +°°)  and 
this  real  function  is  such  that: 
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7lv>(rf 

r  =  [!Z.>.  ’\-da  <  +oo,  Eq.  6.1. 

t  ® 

then 

\\Wf(u,sfdu%. 

—  o~°  s 

This  proof  demonstrates  the  possibility  of  wavelet  decomposition  of  signals  and  is 
predicated  on  the  conditions  described  by  Eq.  6.1,  which  is  called  the  admissibility 
condition.  For  the  integral  described  by  the  admissibility  condition  to  be  finite,  the  wavelets 
must  have  an  average  of  zero. 

6.1  Mexican  hat  wavelet. 

The  wavelets  equal  to  the  second  derivative  of  the  Gaussian  are  Mexican  hat 
wavelets.  The  normalized  Mexican  hat  is: 

2  ft2  > 

m  =  ■ 1  fH"  2^, ' 


Fig.  6.  3  shows  an  analysis  of  the  radar  return  signal  example  (Fig.  6.1).  It  can  be  seen 
that  use  of  this  wavelet  picks  out  the  relevant  features  of  the  signal. 


A.  B. 
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Mocicja  Jut  vavelct:  Tgt:i ,sig:#occi  ,*x:-2t. 9899997, &l:i3. 9099999* 


C. 

Fig.  6.3.  Analysis  of  Radar  example  with  Mexican  hat  wavelet 


6.2  Local  Regularity  of  signals  and  Gaussian  differential  wavelets. 

Local  signal  regularity  is  characterized  by  the  decay  of  the  wavelet  transform  across 
scales.  The  local  regularity  of  a  signal,  f(t),  is  characterized  by  the  Lipschitz  exponents  at  a 
point.  Lipschitz  regularity  for  a  function,  fit),  is  defined:  A  function  /  is  pointwise 

Lipschitz  a  >  0  at  v  if  there  exists  K>  0,  and  a  polynomial  pv  of  degree  m  =  [a  J such  that 
Vr  €  SR,  |  f(t) -  pv{t) |  <  K\t -  v\a  .  Eq.  6.2. 

A  wavelet  transform  estimates  the  exponent  a  if  that  wavelet  has  n  >  a  vanishing  moments: 


J  tkys(t)dt  =0  for  0  <k<n. 


A  wavelet  with  n  vanishing  moments  is  orthogonal  to  polynomials  of  degree  n  -  1 .  Since  a 
<  n,  the  polynomial  pv  has  degree  at  most  n  - 1.  Therefore: 

WpMs)  =  =  0 


and  since  /  =  pv+£v. 


Wf(u,s)  =  Wev(u,s). 

It  also  follows  that  a  wavelet  with  n  vanishing  moments  can  be  written  as  the  /Ith  order 
derivative  of  a  function  6.  Fig.  6.4  below  shows  a  wavelet  transform  of  the  radar  signal  of 

Fig.  6.1  calculated  with  y/  =  -0 ,  where  6  is  a  Gaussian.  The  resulting  Wfiu,s)  is  the 
derivative  of/averaged  in  the  neighborhood  of  u  with  a  Gaussian  kernel  dilated  by  s. 
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It  is  instructive  to  continue  the  analysis  beyond  the  first  derivative,  due  to  the 


following  theorem:  Let  y/  =  (-l)B0(n)  with  0(f)  =  XExp 


r 

2/J2 


•  f  =  fo*8a  and f0  is 


uniformly  Lipshitz  a  on  [v -  h,v  +  h],  then  there  exists  A  such  that: 


V(u,s)e  [v-/i,v  +  /i]x9r,  \Wf(u,s)\  <  Asa+1/2 
where  ga  is  a  Gaussian  variance  of  o2: 


f 


_2  \ 


n-a 


1+- 


pv 


g^,)=^EXP 

Fig.  6.5  below  shows  a  wavelet  transform  of  the  radar  signal  of  Fig.  6.1  calculated 

with  y/  =  -0",  where  0  is  a  Gaussian.  Essentially,  the  above  theorem  describes  how  the 
wavelet  transform  decay  relates  to  the  amount  of  diffusion  of  a  singularity.  The  difference 
between  the  analyses  offered  by  Fig.  6.4  and  6.5  reflects  the  differences  in  the  smoothness 
of  the  decay  of  transitions  in  the  signal.  These  differences  reflect  variations  in  the  signal 
which  can  be  infinitely  continuously  differentiable. 


Gaussian  Differential;  Tgt:i 


0  25  50  75  100  125  150  17S 


A. 


Gaussian  Differential;  Tgt:l ,sig:#00Ol 0899997, Alti3.#099S99* 
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CJUStJan  Dif  f  wrtntlil  r  Tgtll  ,Cigi*»0<U,At!-2«.9(9gs*7,,tril3. 1099994* 


c. 

Fig.  6.4.  Analysis  of  Radar  example  with  Gaussian  differential  wavelet  l ff  =  —0  . 


Gaussian  2nd  Differential:  Tgt:i ,Sig?#oooi ,fci:-2«.«892#$7, *1:13. 8099339* 

r-> - -  - - - - - -  ,  ,  ,  aeusaian  2nd  Differential:  T7t:l  ,$ig:*000l,At:-2«.*e»9i*7,Alii3.80*S2m 


c. 


Fig.  6.5.  Analysis  of  Radar  example  with  Gaussian  2nd  differential  wavelet  If/  =  -6". 

6.3  Spectrogram 

The  windowed  Fourier  transform  requires  a  real  and  symmetric  window  g(t) 
g(-t),  which  is  translated  by  u  and  modulated  at  the  frequency  t,: 


gu4(t)  =  Exp[ify]g(t-u). 

It  can  be  normalized  \\g\\  =  1,  so  that  II gu$  =  1  for  any  (u,%)  e  SR2.  The  resulting  windowed 
Fourier  transform  /  e  L2(SR)  is: 


+0O 

Sf  (U,£)  =  (f,g^)  =  J  f(t)g(t  -  u)Exp[-fy}it, 


also  known  as  the  short  time  Fourier  Transform. 

An  energy  density,  called  a  spectrogram,  Ps,  can  be  defined: 


psfM  =  \sf(u4  = 


J  f(t)g(t  -  u)Exp[-i%t]di 


Fig.  6.6  is  the  spectrogram  of  the  linear  chirp  shown  in  Fig  6.2.  Fig.  6.7  is  the 
spectrogram  of  the  radar  signal  example  of  Fig.  6.1.  It  may  be  seen  that  the  spectrogram 
can  pick  out  well  the  onset  and  offset  of  the  various  modulations  in  the  signal.  However, 
due  to  the  fact  that  the  window  is  constant  (even  although  the  modulated  frequency  of  the 
kernel  is  changing),  the  spectrogram  provides  a  biased  representation  of  the  signal  as  the 
sampling  window  excludes  the  influence  of  signal  components  outside  the  unchanging 
windowed  “box”.  The  spectrogram  is  thus  not  an  optimum  representation  of  signals  and 
future  work  will  demonstrate  how  to  amend  these  methods  to  obtain  wavelet 
representations  from  a  similar  commencement 
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#p*etrc<3r«»:  Linear  CHirp 


c. 

Fig.  6.6.  Spectrogram  or  energy  density  of  a  linear  chirp. 


sp*^10?****-  *gt:i  ,sigi#o*oi  ,fti:-34.»<m»?.*l:i9.«o»»SB9i 


Fig.  6.7.  Spectrogram  or  energy  density  of  radar  signal  example. 


7.0  Zak  Transform. 

The  action  of  the  Zak  Transform  was  investigated  using  a  radar  return  signal 
exemplar  with  variance  of  the  Zak  summing  variable  k.  It  was  shown  that  discrimination 
capability  of  fine  structure  is  lost  as  the  range  of  k  increases.  In  order  to  investigate 
improvements  in  discrimination  capability,  the  Zak  transform  was  performed  not  on  the 
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radar  exemplar  but  on  the  radar  exemplar  after  Gabor  analysis  with  kernels  of  various 
sizes.  It  was  shown  that  there  are  distinctive  patterns  for  each  of  the  kernels  used. 
Therefore  a  radar  return  signal  can  be  analyzed  into  a  4-dimensional  pattern  signature  when 
the  results  of  the  analysis  patterns  of  multiple  kernels  are  summed.  This  method  holds 
promise  for  increasing  fine  discrimination  capability. 

7.1  Zak  Transform  Analysis 

The  Zak  transform1  of  a  function  /  €  L2(9 l)  is: 


Zf(t,co )  =  F{t,co)  =  '£lf(t  +  k)Exp[2mkoo],  (t,co)  e  x  §i 

keZ 

and  satisfies  the  quasiperiodicity  relations: 

Zf(t  + 1,(0)  =  Exp[-2mco]Zf(t,o))  and  Zf(t,co  + 1)  =  Zf(t,co). 

Therefore  the  values  of  the  Zak  transform  are  represented  on  the  cube: 

<2  =  [0,1)  x  [0,1), 

and  the  transform  is  a  unitary  mapping  of  L2(SK)  onto  where: 

L\Q)  =  {f  :  IlFt  =  <  ~j. 

The  action  of  the  Zak  transform  on  a  Gabor  system  {gm  „}  constructed  with  a  =  b  = 
1  is: 

z(8m,J(t,G))  =  Exp[2mmt]Exp[2mna)]Zg(t,co) , 

but  Zg(t,co )  =  1  for  (t,co)  e  Q.  Therefore  the  Zak  transform  maps  the  orthonoimal  basis 
{&n,n}  for  L2( 9i)  onto  the  orthonormal  basis  {Exp[2mmt]Exp[2mnco]}  for  L2(Q). 

1.2  Relation  of  Zak  Transform  to  Gabor  Transform 

The  signal,  f(t),  is  represented  in  a  Gabor  expansion  as: 


1  Daubechies,  I.,  Grossman,  A.  &  Meyer,  Y.  Painless  nonorthogonal  expansions.  J.  Math  Phys  27 
1271-1283, 1986;  ’ 

Daubechies,  I.,  The  wavelet  transform,  time-frequency  localaization  and  signal  analysis.  IEEE  Trans.  Info. 
Theory ,  36, 961-1005, 1990; 

Heil,  C.  &  Walnut,  D.,  Continuous  and  discrete  wavelet  transforms.  SIAM  Review,  31, 628-666,  1989; 
Janssen,  AJ.E.M.,  Bargmann  transform,  Zak  transform,  and  coherent  states.  J.  Math  Phys  23  720-731 
1982; 

Janssen,  A.J.E.M.,  The  Zak  transform:  A  signal  transform  for  sampled  time-continuous  signals  Philips  J 
Res.,  43,  23-69,  1988; 

Zak,  J.,  Lattice  operators  in  crystals  for  Bravais  and  reciprocal  vectors.  Phys.  Rev.  B,  12,  3023-3026, 1975. 
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f(t)  =  ~  mT)Exp[i2jtkUt], 

m  k 


where  the  variable  m  is  used  in  connection  with  a  time  shift,  the  variable  k  in  connection 
with  a  frequency  shift,  and  g(t)  is  the  envelope  of  the  elementary  signal  -  and  conceived  in 
this  case  as  the  synthesis  window 2.  Thus  the  Gabor  expansion  is  the  representation  of  a 
signal  as  the  superposition  of  elementary  packet  signals,  and  the  Gabor  coefficients  a : . 
represent  the  complex  amplitude  of  a  packet  at  time  position  mT  and  frequency  kU. 

With  an  analysis  window ,  y(t),  corresponding  to  the  synthesis  window,  g(t),  the 
Gabor  coefficients,  a^,  are  determined  by  means  of  the  Gabor  transform: 

amk  ~  J  -  mT)Exp[-i2nkUtflt . 

Alternatively,  the  Gabor  transform  can  be  considered  as  a  sampled  version  of  the 
windowed  Fourier  transform  Wft,u)  of  the  signal  f(t): 

Wt(t,u)  =  jf(t  )y*(f  -t)Exp[-i2mf]df , 

on  the  rectangular  lattice  (r  =  mT,  u  =  kU):  a^  =  WfmT,kU). 

There  are  three  sampling  cases  of  note: 

(1)  Critical  Sampling  -  UT  =  1  -  for  which  there  exists  a  unique  relationship  between  the 

synthesis  window,  g(t),  and  the  analysis  window,  y(t),  and  a  synthesis  window  can  be 
uniquely  determined  from  a  given  analysis  window.  In  the  case  of  critical  sampling,  the 
Gabor  signal  expansion  is  related  to  the  degrees  of  freedom  of  the  signal.  If  a  signal  f(t)  is 
limited  to  the  time  interval  |f|  <  0.5  a  and  its  Fourier  transform  to  the  frequency  interval 
\u\  <0.5  b,  the  number  of  degrees  of  freedom  equals  the  number  of  Gabor  coefficients  in 
the  time-frequency  area  ab,  which  is  equal  to  the  time-bandwidth  product  ab.  In  the  case  of 
critical  sampling,  the  representation  is  not  always  stable.  This  instability  is  expressed 
formally  by  the  Balian  Low  Theorem3.  The  stability  problem  can  be  overcome  with 
multiwindows.  A  Gabor  system  is  not  only  a  frame  (see  below),  but  is  also  an  exact  frame 
(see  below)  if  UT  =  1. 

(2)  Oversampling  -UT  <  1  -  for  which  such  a  unique  relationship  does  not  exist  and  there 
is  no  unique  relationship  between  the  synthesis  window  and  the  analysis  window.  A  Gabor 
system  is  a  frame,  but  not  an  exact  frame  when  UT<\. 

(3)  Undersampling  -  UT>  1  -  for  which  the  sequence  of  representation  functions  is 
incomplete  and  therefore  does  not  constitute  a  frame3.  A  Gabor  system  is  incomplete  and 
not  a  frame  when  UT>  1. 

The  Fourier  transform  of  the  expansion  coefficients  is: 

£(&*?)  =  (Fa)(£,7i)  =  X£am*E*pR27r(m77 -/:£)], 

m  k 


2  Bastiaans,  MJ.,  Gabor’s  expansion  of  a  signal  into  Gaussian  elementary  signals.  Proc  IEEE  68  538- 
539,1980. 

3  Daubechies,  I.,  The  wavelet  transform,  time-frequency  localization  and  signal  analysis.  IEEE  Trans  Info 
Theory,  36,  961-1005,  1990. 
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and  it  can  be  shown  that  in  the  case  of  integer  oversampling,  i.e.,  1IU  =  pT,  that4: 


=  pT  +  kpT)Exp[-i2nk(pT)(ripU)} 

k 


X 


X  y{&T  +  [kP~  m]T )Exp[-iljt(kp  -  m)T{r\pU] 

m 


Defining  the  Zak  transform  as: 


Z/I.WA  )  =  £/(»  +  mA)Exp[-i27CmAu], 

m 


then  the  previous  equation  can  be  expressed  using  the  Zak  transform  twice,  with 
t  =  &T ,  u  =  T\pU ,  A  =  pT  for  the /(f)  expression  and  t  =  £pT,  u  =  TjpU,  A  =  T  for 

the  y(t)  expression.  In  terms  of  this  Zak  transform,  the  Fourier  transform  on  the  expansion 
coefficients  is: 


a&ri)  =  pTZf  (&T,  ppU;  pT)Z*r  (&T,  ppU;T). 


7.3  Balian-Low  Theorem 

We  use  the  following  definitions: 

||/||2  =  (J |/(f)|2 dt j  for  the  L2  -norm  of  a  complex- valued  function/,  or  the  energy  of  f, 

L2  being  the  space  of  all  functions  with  finite  energy. 

{/,#)  -  \f(t)g*(t)dt  is  an  inner  product  on  L2(9i). 

A  sequence  {/.}  of  functions  in  L2(9i)  is  orthonormal  if: 


( f&’fi )  —  &k,i 


1,  if  k  =  l 
0,  if  k*l 


A  sequence  {/.}  is  a  basis  for  Z?(9t)  if 


V/  e  L2(5R),  3  unique  scalars  ck(f)  such  that  /  =  Xc*(/)/*>  where  the  basis  elements 

k 

are:  )  =  {f’fk)-  For  an  orthonormal  basis: 


V/  €  /  =  £{/,/,)/, 

Plancherel’s  formula  gives  that  the  energy  of/is  related  to  the  energy  of  the  coefficients: 


4  Bastiaans,  M.J.,  Gabor’s  signal  expansion  in  optics,  pp.  427451  in  Feichinger,  H.G.  &  Strohmer,  T., 
(Ed.s)  Gabor  Analysis  and  Algorithms,  Birkhauser,  Boston,  1998. 
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V/SL2(3i),  |/E  =  £K/./.f 

k 

A  sequence  {/*}  is  &  frame  if  there  exists  AJS  >  0  (the  frame  bounds)  so  that: 


V/Ei?(3i),  A||/g<X|{/,/,>|2£S||/t 

k 

Therefore  for  a  frame  the  energy  ||/|[2  of  f  is  equivalent  or  related  to  the  coefficient  energy 
X|(/,/*>|  .If  A  =B,  then  the  frame  is  right. 

k 

For  any  frame  {/}  there  is  a  dual  frame  j/t  j  so  that: 

v/sm),  /=!(/,/,}/, =£</, /,)/,. 

k  k 

If  the  frame  is  a  basis  then  ck  =  (f,fk^  minimizes  the  energy  ^|c4j2. 

K; 

A  frame  can  be,  but  need  not  be,  a  basis.  If  a  frame  is  a  basis,  then  it  is  exact.  A 
frame  {/}  is  exact  if  there  is  no  single  element  /  which  is  deleted.  A  frame  is  a  basis,  if 

and  only  if,  it  is  exact.  Therefore  an  exact  frame  satisfies  both  the  frame  bound  formula, 
above,  as  well  as  uniquely  representing  a  function /in  terms  of  frame  elements,  fk,  or  dual 
frame  elements,  fk. 

A  Gabor  system  is  defined  by: 


gmb.ruait)  =  Exp[2mmbt\g{t  -  na) . 

With  these  definitions,  the  Balian  Low  Theorem  (BLT )  is: 

Let  {gmb.n*  =  Exp[2mmbt]g(t-na)}m  n£Z  with  a, b  >  and  ab  =  1  (i.e.,  it  is  an  exact  frame). 
As  this  Gabor  system  is  an  exact  frame  for  L2(9i) ,  then 


f  4o°  \/  +oo 

J|rg(r)|2dr  J|)g(r)|2dy 
A— 


=  +oo 


or 


M2lb£(r)||2  =  +~ 


Benedetto  et  al5  point  out  that  the  BLT  maximizes  die  classical  uncertainty  principle,  which 
is: 


Let  (f0>7o)  e  91  x  then 


5  Benedetto,  JJ.,  Heil,  C.  &  Walnut,  D.F.,  Gabor  systems  and  the  Balian-Low  theorem,  pp.  85-122  in 
Feichtinger,  H.G.  and  Strohmer,  T.  (Ed.s)  Gabor  Analysis  and  Algorithms:  Theory  and  Applications, 
Birkhauser,  Boston,  1998. 
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V/  e  L2(SK),  ||/|g  <  4n\\(t  -  ?0)/(0||2||(y  -  7o)/(7)||2 


The  crucial  observation5  is  that  if  a  Gabor  system  {g^  ^}  forms  an  exact  frame,  the  right 

side  of  inequality  is  infinite  when /is  replaced  by  g.  Essentially,  the  BLT  states  that  an 
exact  frame  is  a  filter  set  of  sampling  windows  with  an  infinite  time-bandwidth  product. 
Conversely,  the  uncertainty  product  states  that  there  is  a  minimum  time-bandwidth  product 
for  sampling  windows.  Therefore  the  exact  frame  Gabor  system  can  be  used  to  sample 
signals  of  any  size  time-bandwidth  product,  but  there  is  a  minimum  time-bandwidth 
product  into  which  any  maximum  time-bandwidth  product  can  be  analyzed.  The  condition 
ab  =  1  can  thus  be  interpreted  as  a  Nyquist  criterion. 

7.4  Exemplars 

The  radar  return  signal  exemplar  is  shown  in  Fig.  7.1  (top),  together  with  Gabor 
analyzing  functions,  gm,  for  n  =  1, 2  and  5.  Fig.s  7.2,  7.3,  7.4  and  7.5  show  the  pattern 
analysis  of  the  radar  exemplar  for  the  Zak  transform  variable  k  =  1  to  5,  1  to  10,  1  to  15 
and  1  to  177.  It  can  be  seen  that  discrimination  capability  of  fine  structure  is  lost  as  the 
range  of  k  increases.  Thus,  there  is  a  loss  of  discriminated  fine  structure  in  the  progression 
from  Fig.  7.2  through  to  Fig.  7.5. 

In  order  to  investigate  improving  discrimination  capability,  the  Zak  transform  was 
performed  not  on  the  radar  exemplar  but  on  the  radar  exemplar  after  Gabor  analysis  with 
kernels  of  various  sizes.  Fig.s  7.6,  7.7  and  7.8  show  the  results.  Fig.  7.6  is  the  Zak 
transform  of  the  radar  exemplar  after  filtering  with  the  smallest  Gabor  analyzing  kernel  -  A 
of  Fig.  7.1;  Fig.  7.7  uses  the  next  largest  -  B  of  Fig.  7.2;  and  Fig.  7.8  uses  the  largest  -  C 
of  Fig.  7.1.  It  can  be  seen  that  there  are  distinctive  patterns  for  each  of  the  kernels. 
Therefore  a  radar  return  signal  -  such  as  the  exemplar  -  can  be  analyzed  into  a  4- 
dimensional  pattern  signature  when  the  results  of  the  analysis  patterns  of  multiple  kernels 
are  summed. 


Tgt:l ,Sig:#0001 ,Az: -24. 98999 97,Al: 13, 80999994 
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Fig.  7.6.  Zak  Transform  of  g(m,n),  n  =  1  (Bottom  right),  k  =  177. 
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8.0  Adaptive  Methods 


We  addressed  new  methods  for  achieving  discrimination  of  both  frequency-dependent 
aspects  of  the  target  and  tune-dependent  scattering  centers.  The  Adaptive  Gabor 
Representation  (AGR)  was  used  to  analyze  a  radar  signal  return  in  two  forms:  the  functions 
were  either  varying  Q  or  constant  Q.  We  found  the  constant  Q  form  appears  optimum 
Frequency  hopping  signals  as  well  as  monocycle  pulses  were  analyzed  with  the  AGR  and 
die  Combined  Wigner-Ville  Distribution  of  the  AGR  calculated.  It  was  shown  that  with  all 
frequency  components  combined  the  frequency  components  of  the  signal  can  be  identified 
but  not  the  pulse.  However  when  the  Combined  Wigner-Ville  representation  of  the  AGR 
is  calculated  using  only  the  high  frequency  components  of  the  latter,  it  was  seen  that  the 
monocycle  can  be  detected.  Thus,  these  methods  seems  entirely  appropriate  for  detecting 
the  frequency-dependent  and  the  time-dependent  (transient)  scattering  centers  of  targets. 


8.1  Adaptive  Gabor  Representation  (AGR)6 

The  adaptive  signal  expansion  is  defined: 

s(t)=^Bphp(t), 

P 

with  the  coefficients  defined  by: 

Bp  —  ( s,hp Y 

The  aim  is  to  find  h^(t)  which  is  most  similar  to  s0(t),  i.e., 

M  =  max|(^  (t)hp  (r))|2 .  Step  1 . 

Residuals  are  defined  by: 

Vi  (*)  =  spW  -  Bphp(t)  •  Step  2. 

If  the  functions,  hp(t),  have  unit  energy: 

Mf=i. 

then  the  residual  energy  is: 

6  Qian,  S.,  Chen,  D.  &  Chen,  K„  Signal  approximation  via  data-adaptive  normalized  Gaussian  functions 
and  its  applications  for  speech  processing.  Proc.  ICASSP-92,  San  Francisco,  CA  March  23-26 
1992,  pp.  141-144.  ’ 

Qian,  S.  &  Chen,  D.,  Signal  representation  using  adaptive  normalized  Gaussian  functions  Signal 
Processing,  36,  No  1, 1-11,  March,  1994. 

Mallat,  S.  &  Zhang,  Z.,  Matching  pursuit  with  time-frequency  dictionaries.  IEEE  Trans  Signal 
Processing  41,  3397-3415, 1993. 
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Step  1  is  repeated  to  find  f\  ( t )  as  a  match  to  Sj  (t),  then  Step  2  is  repeated  and  so  on  as  the 
residual  signal  converges  to  zero.  It  should  be  noted  that,  unlike  the  regular  ’  Gabor 

expansion,  as  well  as  wavelets,  the  set  [hp{t)}  will  never  be  complete  in  L2>  even  if  the 

residual  converges  to  zero,  because  each  set  of  adaptive  elementary  functions  is  uniaue  to  a 
particular  signal.  H 

As 

IMf  =  j>,|! 

P= o 


and  as 


±HwVD,r(t,co)dtdm  =  lhp(,f  =  1, 


the  Adaptive  Spectrogram  (AS)  can  be  defined  as: 

AS(t,a»  =  ^\Bp\2WVDhp(t,co), 

P 

where  it  should  be  noted  that  the  summation  is  over  rows. 

The  AS  not  only  does  not  contain  the  cross-term  interference  of  the  Wiener- Ville 
distribution,  it  also  satisfies  the  energy  relation:  & 


k(Of  =  ^-\\AS{t,(o)dtda). 


2  n 

The  choice  of  h/t)  is  completely  arbitrary.  Following  convention?,  we  choose  the 
Gabor- type  functions: 


f  aS^'A 


exp 

\n)  ^ 


where 


so  that: 


At  = 


Af.At  = 


1 


7  Qian,  S.  &  Chen,  D.,  Joint  Time-Frequency  Analysis,  Prentice-Hall,  New  York,  1996. 
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The  Adaptive  Gabor  Representation  (AGR)  is  then: 


s(0  =  X  #/,«>  =  £  S, 

P  P 


An  important  property  of  the  AGR  is  that,  unlike  the  Gabor  expansion,  for  which  the 
analysis  and  synthesis  functions  are  not  identical,  the  AGR  has  the  same  analysis  and 
synthesis  functions.  The  adaptive  coefficients,  Bp,  are  computed  as  follows: 


Bp=j sp(t)hp(t)dt  =  J Sp(t) exp  -  t0p )2  exp[-t/0  r]di 


The  AS  can  then  be  defined  as: 


AS(t,(o)  =  2£ 


where  again  it  should  be  noted  that  this  is  die  summation  of  the  Wigner-Ville  distribution 
over  rows. 

8.3  Wigner-Ville  of  the  AGR. 


The  Wigner-Ville  joint  time-frequency  density  function  of  the  adaptive  Gabor 
functions  is: 


WVDh  (t,  a)  =  2  exp 


-<ccpit-t0f  + 


(CO-Q)0)2 


Testing  the  above  definitions,  the  AGR  was  calculated  of  the  representative  radar 
return  signal  shown  in  Fig.  8.1.  Two  protocols  were  tried  (Fig.s  8.2  and  8.3).  In  the  first 
(Fig.  8.2)  the  AGR  was  calculated  using  5  frequencies,  5  modulating  envelope  widths, 
i.e.,  it  is  not  a  wavelet  representation.  The  center  frequency  is  set,  then  5  modulating 
envelopes  are  applied  in  order  to  ascertain  the  maximum.  In  the  second  (Fig.  8.3)  the  AGR 
was  calculated  using  a  wavelet  representation  in  which  changing  the  center  frequency 
changes  the  bandwidth,  both  temporal  and  frequency.  Then  the  succeeding  analysis  is 
based  on  the  residual.  There  is  indication  that  the  Gabor  wavelet  (Fig.  8.3)  appears  to  be 
the  optimum. 


The  frequency  hopping  signal  (Fig.  8.4)  was  analyzed  using  the  AGR.  Fig.  8.5 
shows  that  the  AGR  can  pick  out  the  individual  signal  frequency  components.  Fig.s  8.6 
and  8.7  shows  the  full  AGR  with  bandwidths  2  x  and  6  x  the  Gabor  wavelet  bandwidth.  At 
6  x  (Fig.  8.7)  there  is  less  precise  timing  identification. 

Next,  a  frequency  hopping  signal  with  a  high  frequency  monocycle  was  analyzed 
(Fig.  8.8).  The  AGR  is  shown  in  Fig.  8.9  and  it  can  be  seen  that  both  the  frequency  and 
the  pulse  components  can  be  identified.  The  Combined  Wigner-Ville  of  the  AGR  of  Fig 
8.9  is  shown  in  Fig.  8.10.  It  can  be  seen  that  the  frequency  components  of  the  Fig.  8.8 
signal  can  be  identified,  but  not  the  pulse.  However,  when  the  Combined  Wigner-Ville 
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representation  of  the  AGR  is  calculated  using  only  the  high  frequency  components  of  the 
latter,  it  is  seen  that  the  monocycle  can  be  detected. 

*gt:l ,sig:#0001  . 9899997,Al:13 . 80999994 


Fig.  8.1.  Representative  radar  return  signal 


s  x  5  *d*ptivC  caber  Kweltt  Tgtn  >«lgt*oooi  ,A»-ai.Basa»B7,AUi3.»09»flsgi  5  X  5  Adaptive  Gabor  xavelet  ygtii.UgitfMOi.i 


ill-21. 989*997  8099 9991 


Fig.  8.2.  Adaptive  Gabor  Representation  (AGR)  of  the  representative  radar  signal  of  Fig.  8.1.  -  5 
frequencies,  5  modulating  envelope  widths.  This  is  not  a  wavelet  representation.  The  center  frequency  is  set, 
then  5  modulating  envelopes  are  applied  in  order  to  ascertain  the  maximum. 


c 

Fig.  8.3.  Adaptive  Gabor  Representation  (AGR)  of  the  representative  radar  return  signal  of  Fig.  8.1.  This  is 
a  wavelet  representation  in  which  changing  the  center  frequency  changes  the  bandwidth  both  temporal  and 
frequency.  Then  the  succeeding  analysis  is  based  on  the  residual.  As  there  seems  little  difference  between 
this  method  and  the  preceding,  the  Gabor  wavelet  appears  to  be  the  optimum. 
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c 

Fig.  8.6.  Adaptive  Gabor  Representation  of  the  frequency  hopping  signal  with  6  x  the  bandwidth  of  the 

Gabor  wavelet.. 
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Adaptive  Gabor  wavelet  frequency  Hopping  Signal*  Pulse  Adaptive  Gabor  Wavelet  Fre<jiencr  Hooping  Signal  +  pulse 


Adaptive  Gabor  Wavelet  Frequency  Hopping  Signal  +  Pulse 


C 

Fig.  8.9  Adaptive  Gabor  Representation  of  the  frequency  hopping  signal  with  monocycle.  Both  the 
frequency  and  the  pulse  components  of  Fig.  8.8  can  be  identified. 
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Combined  Higner-Ville  ranction,  PulseFreqKop  No  i 


Combined  Wigner-ville  ranction,  PulseFreqptop  No  l 


Combined  Kigner-V ille  Function,  PnlseTrecplop,  Ko  1 


Fig.  8.10  The  combined  Wigner-Ville  of  the  AGR  of  the  frequency  hop  signal  with  monocycle.  It  can  be 
seen  that  the  frequency  components  of  the  Fig.  8.8  signal  can  be  identified,  but  not  the  pulse. 
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Fig.  8.11  The  Combined  Wigner-Ville  representation  of  the  Adaptive  Gabor  representation  of  the  frequency 
hopping  plus  monocycle  signal,  using  the  high  frequency  components  of  the  latter.  It  is  seen  that  the 

monocycle  pulse  can  be  detected. 

9.0  The  Constant  Q  Gabor  Wavelet 

The  Adaptive  Gabor  Representation  (AGR)  and  the  Combined  and  Single  Wigner-Ville 
Distribution  of  the  AGR  were  calculated  for  frequency  hopping  signals,  pulses  and  radar 
return  signals  The  AGR  requires  finding  the  optimum  wavelet  center  frequency,  optimum 
dwell  and  optimum  modulating  envelope.  The  hypothesis  under  test  was  that  whereas  the 

AGR  has  been  viewed  as  an  optimization  with  respect  to  three  variables:  ap,  t0  and  f0,  with 

aP  setting  the  envelope  size  independently  of  the  other  two  variables,  use  of  a  constant  Q 

Gabor  wavelet  would  permit  an  efficient  optimization  using  only  one  variable  -  a  -  with 
the  other  signal  parameters  functions  of  that  one  variable.  It  is  shown  that  the  constant  Q 
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Gabor  wavelet  with  f0.t0  =Af.At  =  1/2,  is  optimum  and  that  only  one  variable  is  required 
for  optimization. 

9.0.1  Adaptive  Gabor  Representation  (AGR)8 

The  adaptive  signal  expansion  is  defined: 

P 

with  the  coefficients  defined  by: 


B,  =  (s,hp). 

The  aim  is  to  find  h^(t)  which  is  most  similar  to  s0(t),  i.e., 

\BJ[  =  max|(yp(r)^(0)|  .  Step  1. 

np 

Residuals  are  defined  by: 

Vi W  =  spW  ~  BpKW  •  SteP  2- 

If  the  functions,  hp(t) ,  have  unit  energy: 

Kwf-i. 

then  the  residual  energy  is: 

+P+ i(0f  =  \\sp(t)f  -\bp\2. 

Step  1  is  repeated  to  find  f\  (t)  as  a  match  to  sl  ( t ) ,  then  Step  2  is  repeated  and  so  on,  as  the 
residual  signal  converges  to  zero.  It  should  be  noted  that,  unlike  die  regular  Gabor 

expansion,  as  well  as  wavelets,  the  set  {^(r)}  will  never  be  complete  in  L2,  even  if  the 

residual  converges  to  zero,  because  each  set  of  adaptive  elementary  functions  is  unique  to  a 
particular  signal. 

As 


8  Qian,  S.,  Chen,  D.  &  Chen,  K.,  Signal  approximation  via  data-adaptive  normalized  Canscian  functions 
and  its  applications  for  speech  processing.  Proc.  ICASSP-92,  San  Francisco,  CA  March  23-26 
1992,  pp.  141-144. 

Qian,  S.  &  Chen,  D.,  Signal  representation  using  adaptive  normalized  Gaussian  functions.  Signal 
Processing,  36,  No  1, 1-11,  March,  1994. 

Mallat,  S.  &  Zhang,  Z.,  Matching  pursuit  with  time-frequency  dictionaries.  IEEE  Trans.  Signal 
Processing  41, 3397-3415, 1993. 
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-^JJ  WVDhp  C t,a»dtda )  =  \hp(t) |f  =  1, 
the  Adaptive  Spectrogram  (AS)  can  be  defined  as: 


AS(t,co)  =  £|  B$WVDhr  (r,G», 

P 

where  it  should  be  noted  that  the  summation  is  over  rows. 

,.  .,  ^  AS  not  only  does  not  contain  the  cross-term  interference  of  the  Wiener- Ville 
distribution,  it  also  satisfies  the  energy  relation:  s  vme 

IKOlf  =  AS(t,co)dtdco . 


Gabor-type  Unctions:  *  COmpletdy  arbitra^'  Following  convention9,  we  choose  the 


where 


so  that: 


A  t  = 


Af.At  =  -. 

2 

The  Adaptive  Gabor  Representation  (AGR)  is  then: 


fn  V/4  r 

^(o = YuBPhP( r> = 5X1 -f  )  exp  - h ) 

P  n  \  It  )  2  '  r> 


An  important  property  of  the  AGR  is  that,  unlike  the  Gabor  expansion,  for  which  the 
analysis  and  synthesis  functions  are  not  identical,  the  AGR  has  the  same  analysis  and 
synthesis  functions.  The  adaptive  coefficients,  Bp,  are  computed  as  follows* 


9  Qian,  S.  &  Chen,  D.,  Joint  Time-Frequency  Analysis,  Prentice-Hall,  New  York,  1996. 


(a  Y/4  n  ,"| 


The  AS  can  then  be  defined  as: 


AS(t,co)  =  2^|z?;,|2exp 

p 


where  again  it  should  be  noted  that  this  is  the  summation  of  the  Wigner-Vffle  distribution 
over  rows. 


9.0.2  Wigner-Ville  of  the  AGR. 

The  Wigner-Ville  joint  time-frequency  density  function  of  the  adaptive  Gabor 
functions  is: 


WVDh '  ( t,(0 )  =  2  exp 

«Pfr-*o)2  +  —  — 

; 

. 

i  J 

i_ 

9.0.3  The  Constant  Q  Gabor  Wavelet. 

In  previous  expositions  of  the  AGR,  the  four  signal  variables  of  center  frequency, 
fo>  niidperiod,  t0,  frequency  bandwidth  Af  and  time  bandwidth.  At,  are  separately  defined 
For  example, 


Af  = 

1  J 

At  = 

IT 

1 

2 

KaP) 

with 

fo  =  mid  Af, 
t0  =  mid  At. 


Obviously,  this  does  not  provide  a  constant  Q  wavelet. 

However,  a  constant  Q  Gabor  wavelet  is  achieved  from  the  following  definitions10: 


10  Barrett,  T.W.,  Structural  information  theory.  J.  Acoust.  Soc.  Ant.,  54,  1092-1098, 1973. 
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fo=J&K, 


1 

1 

(32 

ap 

II. 

«, 

V2 

1  ’ 

n 

T 

2  cr 


The  constant  Q  Gabor  wavelet  thus  conforms  to  the  following  condition: 


f0.t0  =  A/ .At  = 


and  the  more  general  condition  is11: 


fo-to  —  A/.Ar  —  —  (2n  + 1)  n  =  0,1,2,. 


9.1  Numerical  Tests. 


The  hypothesis  under  test  was  that  whereas  the  AGR  has  been  viewed  as  an 
optimization  with  respect  to  three  variables:  ap,  t0  an df0,  with  ap  setting  the  envelope  size 
independently  of  the  other  two  variables,  use  of  a  constant  Q  Gabor  wavelet  would  permit 

an  efficient  optimization  using  only  one  variable  -  a  -  with  the  other  signal  parameters 
functions  of  that  one  variable.  6  F 


Using  the  frequency  hopping  signal  plus  a  pulse  shown  in  Fig.  9. 1 ,  a  non-constant 
Q  analysis  was  calculated  with  4/0  s  and  1  original  signal  plus  4  residuals  -  Fig.  9  2  The 
Combined  Wigner-Ville  distribution  picks  out  the  frequency  components  of  die  signal  - 


11  Barrett,  T.W.,  On  vibrating  strings  and  information  theory.  J.  Sound  &  Vibration ,  20, 407-412, 1972a 
- ,  Conservation  of  information.  Acustica,  27,  44-47,  1972b. 

The  defoution  precedence  of  signal  parameters:  sequential  versus  simultaneous  information 
Acustica,  27, 90-93, 1972c. 

..  The  conceptual  basis  of  two  information  theories  -  a  reply  to  some  criticisms.  J  Sound  & 


Vibration,  25,  638-642,  1972d. 


- ,  Analytical  information  theory.  Acustica,  29, 65-67, 1973. 

— 7 — r’  Nonlinear  analysis  and  structural  information  theory:  a  comparison  of  mathematical  and 
physical  derivations.  Acustica,  33, 149-165,  1975. 

- ,  Chi  linearizing  nonlinear  systems.  J.  Sound  &  Vibration,  39,  265-268, 1975. 

_r_,Lmeamy  in  secular  systems:  four  parameter  superposition.  J.  Sound  &  Vibration,  41, 259- 


- ,  Information  measurement  I.  On  Maximum  entropy  conditions  applied  to  elementary 

signals.  Acustica,  35,  80-85,  1975.  y 

- ,  Information  measurement  II.  On  minimum  conditions  of  energy  order  applied  to 

elementary  signals.  Acustica,  36, 282-286, 1975. 

- ,  Structural  information  theory  of  sound.  Acustica,  36, 271-281, 1976. 
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Fig.  9.3;  and  the  individual  Wigner-Ville  distributions  pick  out  the  individual  signal 
frequency  components  -  Fig.s  9.4-9.7. 

Calculating  the  AGR  of  the  frequency  components  of  the  signal  but  with  increasing 
modulating  envelope,  Fig.  9.8,  revealed  that  the  constant  Q  Gabor  wavelet  provided  more 
detail  than  did  longer  envelopes  and  did  not  overlap  the  temporal  length  of  the  frequency 
components  in  the  signal.  This  result  established  that  a  constant  Q  Gabor  wavelets  are 
optimum  filters. 

Using  high  frequency  Gabor  wavelets,  the  pulse  was  extracted  from  the  signal 
shown  in  Fig.  9.1  -  Fig.  9.9;  and  using  the  Combined  Wigner-Ville  distribution 
represented  in  Fig.  9. 10. 


Similar  calculations  were  then  performed  on  the  representative  radar  return  signal  - 
Fig.  9.11.  The  AGR  was  calculated  -  Fig.  9.12  -  and  the  Combined  Wigner-Ville 
distribution  is  shown  in  Figs.  9.13(1)  and  9.13(11)  at  two  levels  of  magnification.  Fig.  9.14 
shows  the  Combined  WV  for  the  high  frequency  components  of  the  signal  and  Figs.  9.15  - 
9.18  show  the  WV  for  decreasing  frequency  components  of  the  signal. 


Fig.  9.1.  Frequency  Hopping  Signal  with  Pulse. 

i  X  \  Adaptive  Gabor  Wavelet,  Frequency  Xqpping  Signal  +  Pulse 
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%  x  4  Adaptive  Gabor  wavelet,  frequency  Mopping  Signal  4  Pnlsi 


Fig.  9.3.  Combined  Wigner-ViUe  Distribution  Transform  of  the  4  x  4  Adaptive  Gabor-Wavelet  Transform 
of  Fig.  9.2. 


#1  Higner-viUe  Function,  t  x  «  Adaptive  Gabor  wavelet.  Frequency 
Mopping  Signal  +  Pulse 
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Fig.  9.4.  Single  Wigner-Ville  Distribution  Transform  of  the  1st  part  of  the  4  x  4  Adaptive  Gabor- Wavelet 
Transform  of  Fig.  9.2. 


#2  vigncr-villc  Function,  l  x  l  Adaptive  Cabor  Wavelet,  Frequency 
Hopping  signal  ♦  Pulse 


#2  wign*r-v lilt  Function, 
Hopping  signal  ♦  Pulse 


»  x  t  adaptive  Gabor  Wavelet,  Frequency 


#2  Higner-ville  Function,  s  x  4  Adaptive  Gabor  wavelet.  Frequency 
Hopping  Signal  ♦  Pulse 


Fig.  9.5.  Single  Wigner-Ville  Distribution  Transform  of  the  2nd  part  of  the  4  x  4  Adaptive  Gabor-Wavelet 
Transform  of  Fig.  9.2. 


#3  Wigner-Villc  run  at  ion,  t  x  4  adaptive  Gabor  Wavelet,  Frequency 
Hopping  Signal  +  Pulsa 

1L - 1  1  1  _  ■  =- 


wigner-vme  r unction,  t  x  t  Adaptive  Gabor  wavelet,  Frequency 
Hopping  Signal  *  Pulae 


wi  Hxgner-V ille  Function,  4x4  Adaptive  Gabor  Vavelet,  Frequency 
Hopping  Signal  +  Pulse 


Fig.  9.7.  Single  Wigner-Ville  Distribution  Transform  of  the  3rd  part  of  the  4  x  4  Adaptive  Gabor-Wavelet 
Transform  of  Fig.  9.2. 

5  Adaptive  Gabor  Xavelet,f  =  1.25,  envelope  1-10,  Frequency  Xcppi  5  *  5  *<*«ptive  Gabor  Wavelet,  f  .  l.so,  envelope  t-io,  rreqnency  Hopping 
igual  -f  Pulse  signal  +  Pulse 


m 


W/a 


.  . .  « 

mm 
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x  5  adaptive  Gabor  Wavelet, f  =  1 .00,  envelope  1-10,  Frequency  Hopping 
Signal  Pulse 


5X5  Adaptive  Gabor  wavelet,!  =  o.?5,  envelope  i-io,  rregwncy  Hopping 
Signal  +  Pulse 


Fig.  9.8.  Adaptive  Gabor  Wavelet  Transform  of  parts  of  signal  of  Fig.  9.1,  with  increasing  modulating 


Fig.  9.9.  Adaptive  Gabor  Wavelet  transform  using  high  frequency  components  of  the  signal  0f  Fig.  9.1 


A 


B 


5X5  Adaptive  Gabor  wavelet  Tgt:l ,sig:#oooi ,Az:-2t. 9899997, A1:13. 80999994 


C 

Fig.  9.12.  Adaptive  Gabor-Wavelet  Transform  of  representative  signal  of  Fig.  9.11.  5  f„s  were  used  and  1 
original  signal  plus  4  residuals. 
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Combined  vigner-ville  7 unction,  3x5  Adaptive  Qabor  Novelet  Tgt:i,3ig: 
#0001  ,Jtz:-2«.  9899997, *1:1 3. 90993934 


(DA.  (I)B. 


combined  Higaer-ville  Function,  5X5  Adaptive  Gabor  Havelet  Tgt:i,sig: 
#0001 ,*X : -24. 9039997, Al:  13. $0999994 


Fig.  9.13(1).  Combined  Wigner-Ville  Distribution  Transform  of  the  5  x  5  Adaptive  Gabor-Wavelet 
Transform  of  Fig.  9.2. 
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c 


Fig.  9.14.  Wigner-ViUe  Distribution  Transform  of  the  5th  level  (highest  fiequency  components) 
Adaptive  Gabor- Wavelet  Transform  of  Fig.  9.2. 


#J  xigner-ViLle  runctioa,  5x3  Adaptive  Gabor  wavelet  Tgt:i  ,sig:#oooi, 
kZ:-2\. 9899997, 11:13.90999994 
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*3  Vi^er-Vill*  Junction,  5  ,x  5  Adaptive  Gabor  Wavelet  *gt:l  rS±m#oooi 
Al: -24. 90SS397,Al:13. 00399334 
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#3  Kigner-viUe  Function,  S  x  s  Adaptive  Gabor  wavelet  Tgtsl  ,sig:#0001 , 
AZ I -2$ . 9899997 ,Al:13.80999994 


Fig.  9. 16.  Wigner-Ville  Distribution  Transform  of  the  3rd  level  (middle  frequencies)  x  5  Adaptive  Gabor- 
Wavelet  Transform  of  Fig.  9.12. 
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02  wilier -rill*  Function,  3x3  Adaptive  C*bor  Wavelet  Tgt:l ,Sigj#oOOl , 
ktl‘2%. >899997, *1:1 3. «09S99»« 


If 


#1  Wigner-Ville  Function,  5x5  Adaptive  Gabor  Wavelet  Tgt:l ,sig:#0001  , 
AZC-24. 98999 97 ,A  1:13. 80999994 


Fig.  9.18.  Wigner-Ville  Distribution  Transform  of  the  1st  level  (lowest  frequency  components)  x  5 
Adaptive  Gabor-Wavelet  Transform  of  Fig.  9.12. 


10.0  Extension  of  the  Gabor  Transform  (CQOW) 

We  initiated  investigations  of  wavelets  based  on  orthogonal  parabolic  cylinder  or 
Weber  functions  of  increasing  level  described  by  the  expansions:  Af.At  =f0.t0  =  1/2  (2 n  + 

1),  n  =  0,1,2, .  (i.e.,  Constant  Q  Orthogonal  Weber  Function  Wavelets  (CQOWs)). 

Using  a  frequency  hopping  and  pulse  test  signals,  these  wavelets  are  able  to  detect  both 
resonance  and  pulse  scattering  features.  Using  scattergram  plots  between  three  levels  of 
wavelets  and  two  scale  dilations,  unique  features  of  both  the  test  signal  and  radar  return 
signals  can  be  demonstrated.  This  research  is  unique  in  two  respects :  using  orthogonal 
wavelets  of  different  levels  and  in  the  scatterogram  cross-comparison  of  three  levels  of 
analysis,  which  can  provide  unique  signatures  for  targets.  In  the  case  of  the  test  signal  of  a 
frequency  hopping  signal,  plus  a  pulse,  it  is  seen  that  discernible  cyclic  patterns  occur 
which  may  be  correlated  with  the  number  of  frequencies  in  the  signal  -  four.  In  the  case  of 
the  radar  retran  test  signal,  it  is  seen  that  recognizable  helical  patterns  occur,  which  may 
provide  a  unique  signature  for  targets. 

In  previous  expositions  of  the  Gabor  transform,  the  four  signal  variables  of  center 
frequency,  f0,  midperiod,  t0,  frequency  bandwidth  Af  and  time  bandwidth,  At,  are 
separately  defined.  For  example, 


A/ 


At  = 
with 

/0  =  mid  Af, 
t0  =  mid  At. 

Obviously,  this  does  not  provide  a  constant  Q  wavelet. 

However,  a  constant  Q  Gabor  wavelet  is  achieved  from  the  following  definitions1: 

fo  = 


h  —  d 


P  9 


-L-L 

'32  a. 


A 

V2  1 


At  = 


11 

'2cr 


The  constant  Q  Gabor  wavelet  thus  conforms  to  the  following  condition: 


1  Barrett,  T.W.,  Structural  information  theory.  J.  Acoust.  Soc.  Am.,  54,  1092-1098,  1973. 
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Uh  =  A/.Ar  = 


A  more  general  condition  has  been  derived2: 

Uh  =  Af.At  =  ~(2n  +  1)  n  =  0,1,2, . 

which  provides  a  range  of  orthogonal  wavelets. 

Specifically,  these  wavelets  are  based  on  modulations  which  are  Weber’s  functions. 
Weber’s  equation  is: 


d2P 

dt2 


+ 


(A-f)Z)  =  0, 


where  %  is  a  dimensionless  independent  variable  and  A,  is  a  dimensionless  eigenvalue.  The 
solutions  to  Weber’s  equation  are  the  parabolic  cylinder  or  Weber-Hermite  functions: 


Dn(t)  =  2~"I2H 


'  1 

1 

1 _ 

wfjexp 

_  4  _ 

where 


Hn(t)  =  (-!)"  exp[r2] 


<9exp[-r2] 

dtn 


are  Hermite  polynomials  satisfying: 


and  the  orthogonality  relationship: 


d2x  .  dx  _ 

— 5-  — 2 1 —  2nx  =  0 
dt2  dt  ~ 


2  Barrett,  T.W.,  On  vibrating  strings  and  information  theory.  J.  Sound  &  Vibration,  20, 407-412, 1972a. 

_ ,  Conservation  of  information.  Acustica,  27, 44-47, 1972b. 

- >  xhe  definition  precedence  of  signal  parameters:  sequential  versus  simultaneous  information. 

Acustica,  27, 90-93, 1972c. 

- ,  The  conceptual  basis  of  two  information  theories  -  a  reply  to  some  criticisms.  J.  Sound  & 

Vibration,  25, 638-642, 1972d. 

_ ,  Analytical  information  theory.  Acustica,  29, 65-67, 1973. 

- ,  Nonlinear  analysis  and  structural  information  theory:  a  comparison  of  mathematical  and 

physical  derivations.  Acustica,  33, 149-165, 1975. 

_ .»  On  linearizing  nonlinear  systems.  J.  Sound  &  Vibration,  39, 265-268, 1975. 

_ ,  Linearity  in  secular  systems:  four  parameter  superposition.  J.  Sound  &  Vibration  41  259- 

261,1975. 

_ ,  Information  measurement  I.  On  Maximum  entropy  conditions  applied  to  elementary 

signals.  Acustica,  35,  80-85, 1975. 

_ ,  Information  measurement  II.  On  minimum  conditions  of  energy  order  applied  to 

elementary  signals.  Acustica,  36, 282-286, 1975. 

_ ,  Structural  information  theory  of  sound.  Acustica,  36, 271-281, 1976. 
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J Hm(t)Hn(t)exp[-t2]dt  =  0  for  m*  n. 


The  Weber  functions,  O’th  to  6th  order,  are: 


Weber  Functions  -  see  Fig.  10.1. 

ORDER  (LEVEL) 

O’th  order,  n- 0 

2  -°nH0 

1  ' 
lV2, 

|exp 

-t2' 

4 

1st  order,  n=l 

2~inHl 

UJ 

exp 

7T 

S3 

2nd  order,  n= 2 

1~2,2H2 

f  t  > 

o/2; 

exp 

'-t2' 

4 

3rd  order,  n=3 

2  ~V2H3 

<  t  \ 

exp 

~-t2  i 

4 

4th  order,  n=4 

2~4/2H4 

(— 
K'fe  ) 

exp 

\-t2~ 
4  _ 

5th  order,  n=5 

2~snHs 

£ 

exp 

’-r2l 

4 

6th  order,  n= 6 

2  -*,2H6 

— 1 
^V2; 

exp 

~-r2l 

4 

The  Weber  function  wavelets,  O’th  to  6th  order,  are: 


Constant  Q  Weber  Function  Wavelets,  see  Fig.  10.2 

ORDER  (LEVEL) 

O’th  order,  n= 0 

2~°I2H0 

mm 

mm 

■ 

exp[i(4/l)r] 

1st  order,  n=l 

2~l,2Hl 

exp 

i 

exp[t'(4/3)r] 

2nd  order,  n=2 

2  ~2/2H2 

exp 

S 

exp[t'(4/ 5)r] 

3rd  order,  n= 3 

2  ~V2H3 

exp 

B 

exp[t'(4/7)r] 

4th  order,  n=A 

exp 

exp[t’(4/9)r] 

5th  order,  n  =5 

2  ~5,2H5 

i 

exp 

g 

exp[t'(4/ll)r] 

6th  order,  n=6 

m 

Si 

exp 

B 

exp[t(4/ll)r] 

At  the  O’th  order,  i.e.,  n  =  0,  the  Weber-Hennite  function  is  the  well-known 
Gaussian  and  the  Weber  Function  wavelet  of  the  O’th  order  is  a  unique  Gabor  function  -  a 
constant  Q  Gabor  function  -  due  to  the  two  uncertainty  restrictions  described  above.  This 
O’th  order  wavelet  is: 


a114 
s(t)  =  -£— 


V2 


exp 


aPi  \2 


1 


exp[i8(r -/„)];  s(f)  =  S7exp[-(2/«,)(/-/0)2]exp[-i(/-/0)/8 


or,  with  a=  ,  . 
p  a 


*(0_^2aeXP~K^r)  exp  ;  S(f)  =  -N/f-exp[-2(^(/  -  /o))2]exp[-ia(/-  /0)/8] 

All  of  the  Constant  Q  Orthogonal  Weber  Function  Wavelets,  which  include  the 
above  as  the  O’th  order  wavelet  are: 


ff 


2~nl2H. 


t~t0 
V  a 

V2 


exp 


rt_-t o 

K_a 


eXP 


i(4/(2(n  +  l)))f  *-k 


a 


,  n  =  0,1, 2,3,. 


69 


Integration  showed  that  these  continuous  wavelets  are  reasonably  well  compactly 
supported. 

10.1  Applications 

Fig.  10.3  shows  the  test  signals.  The  first  is  a  frequency  hopping  signal  with  a 
pulse;  the  second  is  a  representative  radar  return.  Fig.s  10.4  and  10.5  show  the  application 
of  the  Constant  Q  Orthogonal  Weber  Function  Wavelets  (CQOW)  to  the  two  test  signals  for 
the  O’th  order  through  the  6th  order  and  at  two  scale  dilations.  It  can  be  seen  that  CQOWs 
provide  excellent  detection  of  both  signal  frequencies  and  pulse  scattering. 

A  new  form  of  analysis  was  tried  using  the  results  of  the  normalized  CQOW 
coefficients.  Scatterograms  were  calculated  (Fig.s  10.6  and  10.7)  for  the  two  test  signals 
by  parametrically  plotting  the  coefficients  for  the  O’th  order,  1st  dilation  or  01,  the  1st 
order,  1st  dilation  or  1 1,  and  the  2nd  order,  1st  dilation  or  21,  i.e.,  01,1 1,21,  then  the  O’th 
order,  2nd  dilation  or  02,  the  1st  order,  2nd  dilation,  or  12,  and  the  2nd  order,  2nd 
dilation,  or  22,  i.e.,  02,12,22.  These  parametric  plots  were  carried  out  among  the  orders 
O’th  through  6th  and  for  2  scale  dilations. 

In  the  case  of  the  frequency  hopping  signal,  plus  a  pulse,  Fig.  10.6,  it  can  be  seen 
that  discernible  cyclic  patterns  occur  which  may  be  correlated  with  the  number  of 
frequencies  in  the  signal  -  four.  In  the  case  of  the  radar  return  test  signal,  Fig.  10.7,  it  can 
be  seen  that  recognizable  helical  patterns  occur,  which  may  provide  a  unique  signature  for 
targets. 


3rd  Order 


6th  Order 


Fig.  10.1.  Weber  functions,  O’th  to  6’th  order. 
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O’th  Order 


3rd  Order 


1st  Order 


2nd  Order 


5th  Order 


Fig.  10.2.  Wavelets  based  on  Weber  functions,  O’th  to  6’th  order. 
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Fig.  10.3.  A.  Frequency  Hopping  plus  Pulse  Test  Signal.  B.  Radar  Return  Test  Signal. 
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Fig.  10.7.  Parametric  plots  of  0th  Order,  1st  Order,  2nd  Order,  1st  Filters  (01,11,21),  2nd  Filters 
(02,12,22),  etc.  ...(42,52,62)  for  the  Test  SignalJRadar  Return  Signal. 


10.2  Paraunitary  Properties 

Using  a  frequency  hopping  and  pulse  test  signals,  CQOWs  wavelets  are  able  to 
detect  both  resonance  and  pulse  scattering  features.  The  characteristics  of  these  wavelets  are 
based  on  HF  filter  principles.  It  is  customary  to  believe  that  a  discrete  time  scaling  (low 
pass)  wavelet  is  first  necessary  to  obtain  the  discrete  time  (high  pass)  wavelet.  Here  we 
show  that  the  expansion  method  permits  an  estimate  of  the  continuous  time  high  pass 
wavelet  from  the  summed  Fourier  spectra  of  the  expansions,  n  =  1,  2,...  Future  work 
will  examine  whether  repeated  analysis  of  the  original  signal  by  higher  order  wavelets  is 
equivalent  to  the  customary  method  of  using  but  one  scaling  function  and  one  wavelet  with 
signal  decimation. 
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Paraunitary  properties  underlie  not  only  the  perfect  reconstruction  capabilities  of 
certain  wavelets  but  the  relation  of  the  low  pass  (scaling)  and  high  pass  (wavelet)  filters. 
These  properties  are  defined  as  follows3.  Some  definitions  are  necessary,  which  are: 

If  H(z)  is  the  z-transform  of  h(n),  then 


H{z)=  |>(n)z- 

n=— 


If  z  =  exp[t  ©],  then  H(z)  is  the  Fourier  transform  of  h(rt). 

H(z)  =  [H(z)]\ 

that  is,  H(z)  is  the  paraconjugate  of  H(z). 

H*(x)  is  the  complex  conjugate  of  H(x)  and  H,(x)  indicates  that  only  the 
coefficients  are  conjugated. 

Hf(x)  denotes  the  transpose-conjugate  of  H(x). 

Using  the  above,  if  H(z)H(z)  =  dl,  for  some  d>  0  and  for  all  z, 

then  the  filter  has  the  paraunitary  property.  Furthermore,  if  E(z)  is  a  filter  analysis  bank, 
and  R(z)  is  a  filter  synthesis  bank,  then  the  perfect  reconstruction  property  is  expressed  by: 

R(z)  =  cz~KE{z). 

Under  the  paraunitary  condition,  the  analysis  filters  are  then  related  as: 

th(n)  =  -c(-l)%(L-n), 

from  which  relation  has  arisen  the  customary  observation  that  the  discrete  time  high  pass 
filter  must  be  obtained  from  the  discrete  time  low  pass  filter  and  that  knowing  the 
continuous  time  low  pass  filter  does  not  permit  the  definition  of  the  high  pass  filter  - 
continuous  or  discrete  time. 

10.3  Applications 

Here  we  obtain  the  high  pass  filter  in  another  way.  The  Fourier  transforms  of  a 
number  of  Weber  functions  are  shown  in  Fig.  10.3.1.  With  the  n  =  0  function  as  the 
scaling  wavelet  (the  first  shown  in  Fig.  10.3.1)  then  the  higher  order  functions  can  be 
viewed  as  approximations  to  the  high  pass  wavelet.  Thus  the  high  pass  wavelet  can  be 
obtained  by  the  inverse  transform  of  the  summation  of  the  Fourier  spectra  of  the  series 
above  n  =  0.  Fig.  10.3.2A  shows  that  summation  for  the  small  number  of  the  series  shown 
in  Fig.  10.3.1.  Fig.  10.3.2B.  shows  the  limiting  extrapolation. 


3  Cf.  Vaidyanathan,  P.P.  Multirate  Systems  and  Filter  Banks,  Prentice  Hall,  1993; 
Vetterli,  M.  &  Kovacevic,  Wavelets  and  Subband  Coding,  Prentice  Hall,  1995. 
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The  inverse  transform  of  Fig.  10.3.2B  provides  the  high  pass  wavelet  which  is 
shown  in  Fig.  3B  (with  the  scaling  function  shown  in  Fig.  10.3.3A  for  cross  comparison). 
The  low  pass  scaling  function  and  the  high  pass  wavelet  were  then  used  to  analyze  die 
frequency  hopping  signal  plus  pulse  shown  in  Fig.  10.3.4.  Figs  10.3.5  and  10.3.6 
provides  the  results.  It  can  be  seen  that  the  low  pass  and  high  pass  features  are  adequately 
displayed.  It  is  straightforward  to  use  this  scaling  function  and  this  wavelet  in  a  multiscale 
analysis,  utilizing  the  customary  decimation  procedures. 

Alternatively,  the  signal  can  be  analyzed,  without  decimation,  by  utilizing  the 
Weber  functions  wavelet  series,  n  =  0, 1, 2, 3,  4,  5,  6,  in  a  Gram-Schmidt-like  expansion. 
Future  work  will  examine  the  advantages  of  this  approach. 


Fig.  10.3.3.  A:  The  scaling  (low  pass)  Weber  function.  B:  the  (high  pass)  wavelet  from  the  inverse  Fourier 
transform  of  the  Fourier  spectra  of  Fig.  2B. 
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C(1) 


C(2) 


C(3) 


Fig.  10.3.5.  Wavelet  transform  of  the  frequency  hopping  signal  plus  pulse  shown  in  Fig.  4,  using  the  low 
pass  Weber  function  wavelet  shown  in  Fig.  10.3.3A.  A(l)-B(l)-C(l),  A(2)-B(2)-C(2),  A(3)-B(3>C(3))  show 
the  real  part,  the  phase  and  modulus  of  the  wavelet  transform. 


C(l)  C(2)  C(3) 


Fig.  10.3.6.  Wavelet  transform  of  the  frequency  hopping  signal  plus  pulse  shown  in  Fig.  10.3.4,  using  the 
high  pass  Weber  function  wavelet  shown  in  Fig.  10.3.3B.  A(1)-B(1>C(1),  A(2)-B(2)-C(2),  A(3)-B(3)-C(3), 
show  the  real  part,  the  phase  and  modulus  of  the  wavelet  transform. 
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11.0  Instantaneous  Phase 

This  section  addresses  instantaneous  phase  information  obtainable  from  application 
of  wavelet  analysis.  We  show  that  a  detailed  thumbprint  of  a  complex  (target  return)  signal 
can  be  obtained  which  can  be  used  as  data  supporting  a  detailed  wavelet  modulus 
description  of  the  same  signal.  We  also  show  a  result  from  calculating  the  differential  phase 
which  also  provides  supporting  data.  It  is  also  shown  that  there  are  optimum  dilations  for 
providing  correlated  measures. 

We  also  calculated  the  probability  distribution  function  for  the  target  return  signal 
processed  with  a  specific  dilation  of  the  wavelet  and  the  mutual  information  for  two 
dilations.  This  form  of  analysis  also  provides  supporting  data  to  the  modulus  thumbprint  of 
the  signal. 

In  summary,  the  phase  information  and  mutual  information  provide  excellent 
supporting  data  to  that  providing  by  the  wavelet  modulus.  In  combination  a  more  detailed 
analysis  of  target  information  is  possible. 


1.  Instantaneous  Phase 

The  Weber  functions  are  defined  as  extensions  of  the  Gabor  transform  with  the  four 
signal  variables  of  center  frequency,  f0,  midperiod,  t0,  frequency  bandwidth  Af  and  time 
bandwidth.  At,  are  separately  defined.  For  example. 


4f  = 

At  = 


(«,) 
U  J 

L  /  .  \ 

1 

1 

2 

Kap) 

with 

fo  =  mid  Af, 


t0  =  mid  At. 


Obviously,  this  does  not  provide  a  constant  Q  wavelet. 

However,  a  constant  Q  Gabor  wavelet  is  achieved  from  the  following  definitions1: 


/o=#V 


4f  = 
At  = 


1  1 


1  Barrett,  T.W.,  Structural  information  theory.  J.  Acoust.  Soc.  Am.,  54,  1092-1098, 1973. 
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The  constant  Q  Gabor  wavelet  thus  conforms  to  the  following  condition: 

-to  =  A/.  .At  =  — , 

A  more  general  condition  has  been  derived2: 


./o-^o  —  4AAf  —  —  (2n  + 1)  «  —  0,1,2, . 

which  provides  a  range  of  orthogonal  wavelets. 

Specifically,  these  wavelets  are  based  on  modulations  which  are  Weber’s  functions. 
Weber’s  equation  is: 


d2P 
dt 2 


+(x-?)d  =  0, 


where  ^  is  a  dimensionless  independent  variable  and  X  is  a  dimensionless  eigenvalue.  The 
solutions  to  Weber’s  equation  are  the  parabolic  cylinder  or  Weber-Hermite  functions: 


where 


Dn(t)  =  2~n/2HI 


//„(?)  =  (-l)"exp[f2] 


<9exp[-r2] 

df 


are  Heimite  polynomials  satisfying: 


2t—_2nx  =  0 
dt 


2  Barrett,  T.W.,  On  vibrating  strings  and  information  theory.  J.  Sound  &  Vibration,  20, 407-412, 1972a. 

_ ,  Conservation  of  information.  Acustica,  27, 44-47, 1972b. 

- >  The  definition  precedence  of  signal  parameters:  sequential  versus  simultaneous  information. 

Acustica,  27, 90-93, 1972c. 

- ,  The  conceptual  basis  of  two  information  theories  -  a  reply  to  some  criticisms.  J.  Sound  & 

Vibration,  25, 638-642, 1972d. 

_ ,  Analytical  information  theory.  Acustica,  29, 65-67, 1973. 

- ,  Nonlinear  analysis  and  structural  information  theory:  a  comparison  of  mathematical  and 

physical  derivations.  Acustica,  33, 149-165,  1975. 

- ,  On  linearizing  nonlinear  systems.  J.  Sound  &  Vibration,  39, 265-268, 1975. 

- ,  Linearity  in  secular  systems:  four  parameter  superposition.  J.  Sound  &  Vibration  41  259- 

261,  1975.  ’  ’ 

- ,  Information  measurement  I.  On  Maximum  entropy  conditions  applied  to  elementary 

signals.  Acustica,  35,  80-85,  1975. 

- ,  Information  measurement  II.  On  minimum  conditions  of  energy  order  applied  to 

elementary  signals.  Acustica,  36, 282-286, 1975. 

_ .  Structural  information  theory  of  sound.  Acustica,  36, 271-281, 1976. 
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and  the  orthogonality  relationship: 


| Hm (t)Hn (r) exp[-r2]dr  =  0  for  m*n. 


The  following  analysis  is  based  on  the  O’th  order  (n  =  0)  wavelet: 


4  J 

The  signal  analyzed  is  the  frequency  hopping  plus  pulse  signal  shown  in  Fig.  1. 
The  wavelet  can  be  viewed  as  an  analytic  signal: 

V r(t)  =  s(t)  +  is(t), 

where  s(t)  is  the  Hilbert  transform  of  sff): 


r(D  =  2^'2H0^jexp 


S(t)  =  -P.vA^-dr, 
n  ij-'t 

P.V.  indicates  that  the  integral  is  taken  in  the  sense  of  the  Cauchy  principal  value.  The 
instantaneous  phase  is  then  defined  as: 


4>(t)  =  arctan— — . 

s(t) 

The  differential  of  instantaneous  phase  is  represented  by  means  of  obtaining  die 
difference  in  phase  between  the  phases  obtained  from  applying  the  Gabor  wavelet  at 
different  dilations. 

11.1  Mutual  Information. 

If  the  phases  of  the  wavelet  filtered  signal  are  considered  as  n  discrete  random 
variables,  Xif...,Xn,  with  sets  of  values,  Ei,...En,  the  probability  distribution  for  an 

individual  phase  for  a  particular  time  in  the  signal  and  using  a  particular  dilation  of  a  Gabor 
wavelet,  Xt,  is: 


P(*i)  =  Pr[*.  =  xt\  xt  €  St.. 

The  joint  probability  distribution  can  then  be  defined  for  n  variables  X,  X  as 

The  mutual  information,  I(Xj;X2),  is  defined  as: 

/(X,;...;X„)  =  X-  Xp(*„~.*.)l°g  '■ 

*,ec,  p{x{)...p(Xn) 
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11.2  Applications 

Fig.  11.2  shows  the  instantaneous  phase  information  from  application  of  the  Oth 
order  Gabor  wavelet  to  the  signal  of  Fig.  1 1.1.  It  can  be  seen  that  a  detailed  thumbprint  of 
the  frequency  hopping  plus  impulse  signal  can  be  obtained  which  could  be  used  as  data 
supporting  a  detailed  wavelet  modulus  description  of  the  signal.  Fig.  11.3  shows  the 
differential  phase  result  which  also  provides  supporting  data.  Fig.  11.4  is  a  more  detailed 
single  cut  analysis  of  Fig.  11.3  and  shows  the  fine  structure  correspondence  to  the 
frequency  components  of  the  signal. 

Fig.  11.5  provides  3-dimensional  cross-phase  plots  of  the  data  of  Fig.  11.2  and 
reveals  the  fine  structural  changes  in  phase  at  increasing  dilations  of  the  wavelet.  Fig.  1 1.6 
also  provides  the  same  information  for  the  differential  phase  measurements  of  Fig.  1 1.3. 
Fig.s  11.7  and  11.8  are  phase-phase  plots  for  different  dilations  of  the  wavelet  for  the 
instantaneous  and  differential  phase  measurements.  It  can  be  seen  that  there  are  optimum 
dilations  for  providing  correlated  measures. 

Fig.  11.9  shows  examples  of  a  probability  distribution  function  for  the  signal 
processed  with  a  specific  dilation  of  the  wavelet  and  the  mutual  information  for  two 
dilations;  and  Fig.  11.10  shows  the  mutual  information  from  application  of  the  Oth  order 
Gabor  wavelet  at  a  number  of  dilations.  This  form  of  analysis  also  provides  supporting 
data  to  the  modulus  thumbprint  of  the  signal. 

In  summary ,  the  phase  information  and  mutual  information  provide  excellent 
supporting  data  to  that  providing  by  the  wavelet  modulus.  In  combination  a  more  detailed 
analysis  of  target  information  is  possible. 


Fig.  11.1.  Frequency  Hopping  plus  Pulse  representative  signal. 
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Fig.  1 1 .3.  Differential  phase  measurements  of  the  signal  of  Fig.  11.1  using  the  Oth  order  Gabor  wavelet 

function. 
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12.0  Fractals  and  Fano  Factors 

This  section  addresses  the  differentiation  of  two  target  return  signals  nsina  m 

Thu^fra6^  ^  iFa^°  faCt°r  pattems  generated  permitting  the  differentiation  of  the  signals 
be  SSSSS^StaK^1 Pr°miSin8  aPPIOaCh  *°  signals  anf wfli 


on 


Q<r> « 


jogQO^I 
log(a)  J 


Q(r)  =  Brf 
V 

where  B,  b  and  a  are  constants  andf(x)  is  a  periodic  function  such  that  f(\+x)  -  fir) 
In  its  simpler  form,  this  relationship  is:  I{  l+x)  -f(x). 


Q(r)  =  Br 

.^S£S£^“S“=:  ffi 

pieces  JS8"8  Aat  a  number  of 


N(r)  »  r~d. 


The  scaling  exponent,  b,  means  that  the  value  of  the  property  O(r) 
r  is  proportional  to  ib\  v  F  h 


measured  at  resolution 


UKr) «  r 


StoS5£?15SfS^2' is  2150  *** ,0  ,he  nnmber' N<r>- of  sa“pks-  ™mpw 


Therefore,  by  substitution: 


and 


Q(r)  =  rN(r). 

Q(r)  =  r.r~d  =  rl~d  =  rb. 
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d  =  l  +  b. 


In  general. 


and  therefore: 


N(r)~r~d, 
Q(r)  ~  rb, 
Q{r)~[N{r )] 


\°  rP 


d  = 


P-b 


a 


where  theoretical  justification  is  needed  for  how  Q(r)  depends  on  N(r>  and  r.  Below  we 
assume  a  =  j8  =  1. 

Another  measure  of  ftactal  behavior,  the  Fano  factor,  is  equal  to  the  variance/m^n 
Sfe^inS6  1S  d6pendent  on  the  window  IenSth  over  which  the  variance  and  mean  are 

•  •  ^  connection  of  ftactal  behavior  and  wavelets  is  due  to  wavelet  coefficients 

giving  local  information  about  differentiability  and  Holder  continuity  of  a  foSn 
Furthermore,  an  essential  characteristic  of  the  wavelet  transform  is  the  Zhilin/  m 
process  relative  to  scale  or  resolution.  Also,  the  wIX  ^s?Su  by  c^ction  tdds 

by  succ!ssive  refinements,  commencing  torn  a  Srs^  aZStiin 
£££&£  ^  SKP  WHh  ^  *““•  ‘  a I™«  identical  TOTtaS 

•  .  Therefore,  we  explored  this  relationship  to  determine  whether  this  form  of  analvsis 

might  provide  a  useful  thumbprint  of  target  return  signals.  ^ 

12.1  Applications 

Fig  miA3115 B^^The8^^0^  °?  tW0  rePresentative  target  signal  returns  shown  in 
rig.  iz.i(A,  B).  The  Gabor  wavelet  transforms  of  these  signals  uDon  which  th* 

subsequent  analysis  was  performed,  are  shown  in  Fig.  12.2.  A  log-log  plot  (amplitude 
124iSS™  scale)  is  shown  for  both  signals  in  Fig.  12.3.  Cuts  throtgluhese  plots  (Fit 
12.4)  show  a  stable  distribution  (as  opposed  to  asymptotically  Gaussian  distribution! 

i^a  fractal  property8  “  ^  ^  °f ' °ne  signa1’ and  decreasing  in  the  case  of  the  other.  This 

was  rr0 

8&E£ XnsZptaf ofKg.Tlt  °f  Rg'  121  overlaid 

tire  signal  ^  *  ****  &ffcrenct  in  Ae  pattemf5 je^S^Sg^e^^n^tiot  of 
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Fig.  12.1.  A:  Representative  radar  return  signal  #1;  B.  Representative  radar  return  signal  #2. 


Fig.  12.2.  (Gabor)  wavelet  transforms  of  signal  #1  (A)  and  signal  #2  (B). 


so  75  100  125  150  175 


Fig.  12.3.  Log-Log  plots  of  wavelet  transforms  shown  in  Fig.  2. 


Fig.12.7.  Fractal  dimensional  representation  of  the  signals  of  Fig.  12.2  with  overlay  of  those  signals. 


Fig.  12.8.  Fano  number  representation  of  the  signals  of  Fig.  12.2. 
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13.0  Group  Theoretical  Approaches 

HitwJSft  86011011  !^esses  in  terms  of  (a)  higher-order  symmetries;  (b) 

JfrS  Ti  (?^ouP  theoretlc’  formulations  in  electromagnetic  theory  of 
die  rde  of  the  A-field  in  (a)  nonlinear  collective  excitations  in  a  2-dimensional  electron  gas* 
and  (b)  sohtons  in  new  SAW  devices  and  heterostructures  in  general.  The  method  of 
generation  of  modulated  TeraHertz  pulses  utilizes  a  plasma  mirror  dependent  on  the 
formation  of  sohton  transmission.  Therefore  the  following  describes  the  commencement  of 
an  analysis  of  such  sohton  transmissions  using  the  various  approaches  described. 

The  equation  of  motion  derived  by  Harmuth1: 


d2E 

dy2 


■He 


d2E 


dE 


dt 


2  -{n°+es)—-saE  =  0 
at 


mcoyorating  5,  the  magnetic  conducmtity,  was  shown  by  Barrett^)  be  a  two-dimensional 
Klein-Gordon  equation  m  the  sine-Gordon  form,  i.e.,  a  soliton  form.  Soliton  solutions 

^^P^mFTle!ei-inti^ablhty-Mdmiegrable  systems  conserve  geometric  features  related  to 
SiUSSc^fTTm6  thC  equatl°nS  of  motion  for  conventional  Maxwell  theory,  which  are 
solutions  of  U(l)  symmetry  systems,  sohtons  are  solutions  of  SU(2)  symmetry  systems 

These  notions  of  group  symmetry  are  more  fundamental  than  differential  aquation 
descnpuons.  Therefore,  this  we  develop  some  basic  concepts  in  order  to  place  differential 
descriptions  within  the  context  of  group  theory. 

Within  this  context,  ordinary  differential  equations  are  viewed  as  vector  fields  on 
mamfolds  or  configuration  spaces.  For  example,  Newton’s  equations  ate  second  “to 
differential  equations  describing  smooth  curves  on  Riemannian  manifolds.  Noether’s 

theorem  states  that  a  diffeomoiphism,  <p,  of  a  Riemannian  manifold,  C,  induces  a 

diffeomorphism,  Dtp,  of  its  tangent  bundle,  TC.  If  0  is  a  symmetry  of  Newton’s  equations, 
then  Dtp  preserves  the  Lagrangian,  i.e., 


Zv°  utp  =  L. 


As  opposed  to  equations  of  motion  in  conventional  Maxwell  theory,  soliton  flows 
are  Hamiltonian  flows.  Such  Hamiltonian  functions  define  symplectic  structures  for  which 


1  Harmuth  H.F.,  Correction  for  Maxwell's  equations  for  signals  I.  IEEE  Trans.  Electromagn  Compat 
EMC-28, 250-258, 1986a;  ’ 

- 1986?WeI1S  eqUati°nS  fOT  SignalS  n' IEEE  TranS-  ElectromaSn-  ComPat., 


Barrett,  T.W.,  Comments  on  the  Harmuth  Ansatz:  use  of  a  magnetic  current  density  in  the  calculation  of 
the  propagation  velocity  of  signals  by  amended  Maxwell  theory.  IEEE  Trans  Electromasn 
Compat.,  30, 419-420,  1988.  S  ' 

- — .  Comments  on  ‘Solutions  of  Maxwell’s  equations  for  general  nonperiodic  waves  in  lossy 

media’”,  IEEE  Trans.  Electromagn.  Compat.,  31, 197-199, 1989. 

- -  Comments  on  ‘Some  comments  on  Harmuth  and  his  critics’”.  IEEE  Trans  Electromasn 

Compat.,  31,  201-202,  1989.  *  ' 

- -  Electromagnetic  phenomena  not  explained  by  Maxwell’s  equations,  pp.  6-86  in  A. 

Lakhtakia  (Ed.)  Essays  on  the  Formal  Aspects  of  Electromagnetic  Theory,  World  Scientific,  1993. 
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“  “d  cha 


_  o  •  ^!i?r<Jer  t0  re*ate  ^e  t^lree  rnaj°r  soliton  equations  to  group  theory  it  is  necessarv  tr> 
examine  the  Lax  equation  or  the  zero-curvature  condition  (ZCC).  ThtZCC  expresses^ 

fpaMs!  ?997)°nnCCtl0n  by  ±e  commutation  relations  of  the  covariant  derivative  operators 


=  0, 


or 


xsThor  that  ^  generic  cases  of  soliton  -  the  Kortewee  de  Vries 
^GE)  cSvene  Norimear  Schrochnger  Equation  (NLS),  the  Sine-Gordon  Equation 
(  °b)  be  ^ven  811  SU(2)  formulation.  In  each  case,  below,  V  is  a  one-dimensional 

space  that  is  embedded  in  the  space  of  off-diagonal  complex  matrices,  f°  ^ 


i  c  0 


and  in  each 


case  A(u)  =  aX  +  u,  where  u  is  a  potential,  A  is  a  complex  parameter,  and  tt  is  the 
constant,  diagonal,  trace  zero  matrix  a  =  f  *  ^ 


SU(2)  formulation. 


0 


which  links  these  equation  to  an 


From  inverse  scattering  theory,  a  function  is  needed,  defined: 


*  1  4~ 

B(®  =  Zc»  exPt-^a  +  —  J b(k)&xp[ik^]dk,  where 
2  2 

are  discrete  eigenvalues  of  u, 

Ci,—CN  are  normalizing  consts,  and 
b(k)  are  reflection  coefficients. 


In  the  first  case,  if  u(x )  = 


B(u)  =  a?i  +  u)i  + 


( o  q(xf 

-l  o  J 

a 

2q 


and 

i  > 
2qx 

ft. 

A  + 

4  2 

i 

~~q 

q 

-qx 

2  J 

u 

4  J 

then  ZCC  is  satisfied  if  and  only  if  q  satisfies  the  KdV  in  the  form  qt 


~t(6 


Ill 


f  o  q(*y\ 

y-qix)  0 


In  the  second  case,  if  u{x)  = 

B(u)  =  atf  +  uX2  + 


and 


then  ZCC  is  satisfied  if  and  only  if  q(x,t)  satisfies  the  NLS  in  the  form 

?«= ■£■(?»+ 2M24 


( 


In  the  third  case,  if  u(x)  = 


B(u)  = 


0 

qx(x) 


4  Ax) 


0 


and 


i 

41 


cos[^(x)]  sin[^(jt)]  \ 
sin[<7(*)]  -cos[#(jt)]J’ 


then  ZCC  is  satisfied  if  and  only  if  q  satisfies  the  SGE  in  the  form  qt  =  sin[$]. 

coin™  Jem  equation  of  moti°n  with  electric  and  magnetic  conductivity  is  in 

Sc  • f0rm’  *e  §r°Up  for  which  is  SU(2).  Solitons  define  Hamiltonian 

flows  and  their  energy  conservation  is  due  to  their  symplectic  structure. 


,  11/v  The  foUowing  are  minimal  topological  and  group  theory  definitions  required  to 
Mow  die  above-initiated  arguments.  The  groups  used  in  Yang-Mills  theory  are  continuous 
groups  (as  opposed  to  discrete  groups).  Unlike  discrete  groups,  continuous  groups  contain 
an  infinite  number  of  elements  and  can  be  differentiable  or  analytic.  V 


U(n)  Group  Algebra 

Unitary  matrices  have  a  determinant  equal  to  ±1.  The  elements  are  represented  bv 
nxn  unitary  matrices.  v  y 


U(l)  Group  Algebra 

Throne-dimensional  unitary  group,  or  U(l),  is  characterized  by  one  continuous 
parameter.  U(l)  is  also  differentiable  and  the  derivative  is  also  an  element  of  U(l).  A  well- 
known  example  of  a  U(l)  group  is  that  of  all  the  possible  phases  of  a  wave  function,  which 
are  angular  coordinates  in  a  2-dimensional  space.  When  interpreted  in  this  way  -  as  the 
internal  phase  of  the  U(l)  group  of  electromagnetism  -  the  U(l) group  is  merely  a  circle  (0 


Mobius  Group,  M,  Algebra 

.  Mobius,  or  M,  transformations  map  the  extended  complex  plane  one-to-one  onto 
itself.  An  element  is 


m:C  ->C  or 

,  \  az  +  b  . , 

m(z)  = - with 

cz  +  d 

ad-bc  *  0. 
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There  is  a  homomorphism  between  M  and  GL(2,C).  The  mapping  from  GL(2,C)  to  M  is: 


p :  X 


a 


b\ 


— »  m(z)  = 


qz  +  ft 
cz  +  d 


The  determinant  of  A 


d 


is  A2det 


,  whence  it  is  seen  that  (1)  1  may  be  chaosen 


in  two  ways  and  that  (2)  the  determinant  can  be  qual  to  +1  -  which  gives  the  grout) 
SL(2,C). 

Because  of  (1)  it  is  seen  that  whereas  SL(2,C)  is  simply  connected,  M  is  not. 
Therefore  SL(2,C)  is  the  universal  covering  group  of  M.  If  u  and  v  are  complex  numbers 
and 


u'  =  au  +  bv 
v'  =  cu  +  dv 


and  if  z  =  u/v  and  w  =  u’lv’,  then 


iv 


az  +  b 
cz  +  d 


a,b,c,d  e  91 


and 


ad-bc  =  1. 


O(n)  Group  Algebra 

The  orthogonal  group,  O(n),  is  the  group  of  transformation  (including  inversion)  in 
an  n-dimensional  Euclidean  space.  The  elements  of  O(n)  are  represented  by  nxn  real 
orthogonal  matrices  with  n(n  - 1)/2  real  parameters  satisfying  AA£  =  1. 

0(3)  Group  Algebra 

The  orthogonal  group, 0(3),  is  the  well-known  and  familiar  group  of 
transformations  (including  inversions)  in  3-dimensional  space  with  3  parameters,  those 
parameters  being  the  rotation  or  Euler  angles  (a,/3,y).  0(3)  leaves  the  distance  squared, 
xz  +  y2  +  z2,  invariant. 


S0(2)  Group  Algebra 

The  collection  of  matrices  in  Euclidean  2-dimensional  space  (the  plane)  which  arc 
orthogonal  and  moreover  for  which  the  determinant  is  +1  is  a  subgroup  of  0(2).  30(2)  is 
the  special  orthogonal  group  in  two  variables. 

The  rotations  in  the  plane  is  represented  by  the  30(2 )  group: 


R(d)  = 


'cos[0] 

vsin[0] 


-sinf#]^ 

cos[0]  J 


0e9? 


where 


R(e)R(y)  =  R(9  +  y). 
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“  ±e  complex  plane  with  multiplication  as  the  group  operation  is  an 
SO(3)  Group  Algebra 

The  collection  of  matrices  in  Euclidean  3-dimensional  snarv»  „„  , i  , 

-dhrrTer  f°r  wh!ch  the  determinant  is  +1  is  a  subgroup  of  0(3).  SO(3)  is  the  special 
orthogonal  group  m  three  variables  and  defines  rotations  in  3-dimensional  space.  ^ 

Rotation  of  the  Riemann  sphere  is  a  rotation  in  9t3  or  £  -  77  -  £  space  for  which 


2x 


N2+i*  77  N2+r 

=*+/y=i±i2 

i-c 


2y 


c= 


z|2-l 

N2+r 


-i\/V°/2 


=^(i  iX 


0  ^  1  r  1  r 


0 


ia  /  2 


tf„03)  = 


1  fl  -A/V*'2  0 


V2 


-/  1 


0 


-iP/2 


rcos  a/2  isin  a/2" 
v* sin  a/2  cos  a/2 ^ 

cos  p/2  -sin  )3/2Y 


£/f(r)  ■  f1  °Y^'2  0  wi 


-J2{-1  1 

\L  f1 

t/2  U  lj  (sin  /i/2  cos  /i/2 

f cos  7/2  -siny/2 
siny/2  cosy/2 


or  ±C/|(a)->/?1(a), 
or  ±Uv(/3) ->  R2(/3), 
or  ±f/f(y)  -»  /?3(y). 


?  rr wi,h 

U^a)  =  ei(an)a', 

Uv(P)  =  e~i(fi,2)a\ 

U((y)  =  eii7,2)ff\ 

(0  n 

cr,  = 

U  0, 

<Ti,<t2,<t3  are  the  Pauli  matrices. 


a2  = 


where 

(0  -r 

0  o. 


<7,  = 


n  o" 
0  -1 


SU(n)  Group  Algebra 

Unitary  matrices  have  a  determinant  equal  to  ±1.  Special  unitarv  matrices  are 
elements  of  unitary  matrices  which  leave  the  determinant  equd  to  +1  ThS  l 
independent  parameters.  SU(n)  is  a  subgroup  of  U(n)  for  whfch  the  detemdS  ^  +1 
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SU(2)  Group  Algebra 

SU(2)  is  a  subgroup  of  SL(2,C).  The  are  22  - 1  =  3  independent  parameters  for  the 
special  unitary  group  SU(2)  of  2x2matrices.  SU(2)  is  a  Lie  algebra  such  that  for  the 
angular  momentum  generators,  JL,  the  commutation  relations  are 

=  tegkJkJJjk  =  L2,3.  The  SU(2)  group  describes  rotation  in  3-dimensional  space 

with  2  parameters  (see  below).  There  is  a  well-known  SU(2)  matrix  relating  the  Euler 
angles  of  0(3)  and  the  complex  parameters  of  SU(2)  is: 


COS 

Ji(cc  +  Y)~ 

'l' 

’-(a-y)' 

.2JexpL  2  J 

sin 

2 

exp 

2 

-sin 

r»r\c 

T 

exp 

'-i(a  +  y)~ 

L2jeXPL  2  . 

.2. 

.  2 

where  a,/3,y  are  the  Euler  angles.  It  is  also  well  known  that  a  homomorphism  exists 
between  0(3)  and  SU(2),  and  the  elements  of  SU(2)  can  be  associated  with  rotations  in 
0(3);  and  SU(2)  is  the  covering  group  of  0(3).  Therefore,  it  is  easy  to  show  that  SU(2) 
can  be  obtained  from  0(3).  These  SU(2)  transformations  define  the  relations  between  the 
Euler  angles  of  group  0(3)  with  the  parameters  of  SU(2).  For  comparison  with  the  above, 

if  the  rotation  matrix  R( a,P,y)  in  0(3)  is  represented  as: 


f  cos  [a]  cos[/3]cos[y]  -  sin[a]sin[y]  sin[a]  cos[/3]  cos[y  ]  +  cos  [a]  sin[y]  -sin[/J]cos[y] 
-cos[a] cos[/3] sin[y]-sin[a]cos[y]  -sin[a]cos[/J]sin[y]  +  cos[a]cos[y]  sin[/J]sin[y] 
v  cos[a]sin[j5]  sin[a]sin[/J]  cos[j3] 

then  the  orthogonal  rotations  about  the  coordinate  axes  are: 


Rl(cc)  = 


( cos[a]  sin[a] 
-sin[a]  cos[a] 


l,  0  0 


0 


f cos[j3]  0  -sin[/3;h 


Rs(r)  = 


0  1  0 

vsin[j3]  0  cos[/3]  j 

r  cos[y]  sin[y]  0^ 
-sin[y]  cos[y]  0 


0  0  1J 


An  isotropic  parameter,  GJ,  can  be  defined: 


GJ  = 


x-iy 

z 
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where  x,y,z  are  the  spatial  coordinates.  If  CJ  is  written  as  the  ..  a 

homogeneous  coordinates  of  the  bilinear  transformation,  then:  Q  M  ^ 


or  the 


|A/4)  = 

cos 

3-Ft2] 
“  />  1  r  .. 

•  \P~ 

sm  — 

L2. 

exp 

1 - 1 

1 

i - 1 

-sin 

PY  f  i(a-y) 

LtH  2 

PI 

cos  £-  exp 

^  Zr 

r-i(a  +  y)  1 

2  J 

mW 


which  is  the  relation  between  the  Euler  angles  r>f  /O/?  i  „ nA  «.v,_  , 

SU(2).  However,  there  is  not  a  unirme  /£,,»  ^  ^  5°?  the _  complex  parameters  of 

correspond  to  1  direction  in  SU(2)  ^here  is  rh?!?  „reanon’  for  2  rotations  in  0(3) 
between  0(3)  and  SU(2).  '  K  S  *  S  a  many‘to'one  or  homomorhpism 


In  the  case  of  a  complex  2-dimensional  vector  (u,v): 


cos 

sin 

T 

exp 

exp 

‘-(a -7)  I 

\ 

-sin 

-2j  *1  2  J 

cos 

2 

r£~ 

.  2 
"-/(a +  7)' 

L2J  L  2 

_2_ 

.  2 

J 

If  we  define: 


then 


a  =  cos 

exp 

~i(a  +  y)~ 

.2. 

2 

b  =  sin 

T 

exp 

'-(a  -7)1’ 

2_ 

-  2  J 

where 


a  b 
-b*  a * 


are  the  well-known  SU(2)  transformation  rules 
determinant: 


Defining:  c  — b'  and  d  =  a,  we  have  the 


ad  -  be  =  1 
or 

aa*-b(-b*)  =  1. 
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(u  v)  as0661™2  (X,:y'Z)  coordinates  with  aspect  to  a  complex  2-dimensioanl  vector 


z  =  uv 

then  SU(2)  transformations  leave  the  squared  distance  x2  +  y2  +  z2  invariant. 
Every  element  of  SU(2)  can  be  written  as: 


Now  if  we  define: 


,  where 


\af  +  \b\2  =  l. 


a  =  yi~iy2> 
b  =  y3  -%» 

then  the  parameters  indicate  positions  in  SU(2)  with  the  constraint: 


y2+y22+y2+yl  =  i. 


E^cates  that  the  group  SU(2)  is  a  3-dimensional  unit  sphere  in  the  4-dimensional  y- 

on  sphere  ***  *>  «■£ 
SU(2)  is  the  quantum  mechanical  “rotation  group”. 


SU(3)  Group  Algebra 

SU(3)  has  n-l  =  8  generators. 

SL(2,C)  Group  Algebra 


SL(2  SPeCial  Unear  of  2  x  2  matrices  of  deteiminant  1  with  comple 


entries  is 


GL(n,R)  Group  Algebra 

The  collection  of  real  non-singular  nxn  invertible  matrices  is  GL(nR)  or  the  real 
general  linear  group  in  n  variables.  (  W  e  1631 

SL(n,R)  Group  Algebra 

SL(n,R)  is  a  subgroup  of  GL( n,R)  with  determinant  =  +1. 


GL(n,C)  Group  Algebra 

The  collection  of  complex  non-singular  nxn  invertible  matrices  is  GUn  C)  or  the 
complex  general  linear  group  in  n  variables. 
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SL(n,C)  Group  Algebra 

SL(n,C)  is  a  subgroup  of  GL(n,C)  with  determinant  =  +1. 

Sp(2,n)  Group  Algebra 

A  2n-dimensional  space  can  be  defined  as: 

Cx,y)  =  (x1y"+1  -  xn+lyl )  +  (x2yn+2  -  xn+zy2)  + ...  +  (xny2n  -  x2nyn) 
-  x'Jy 


Sp(2,n)  is  a  noncompact  symplectic  group. with  form  invariance: 


J  = 


where  1  n  =nxn  unit  matrix. 


Sp(2,n)  satisfy: 


A‘JA  =  J. 


Lie  Group  Algebras 

If  a  topological  group  is  a  group  and  also  a  topological  space  in  which  group 
operations  are  continuous,  then  Lie  groups  are  topological  groups  which  are  also  analytic 
manifolds  on  which  the  group  operations  are  analytic. 

In  the  case  of  Lie  algebras,  the  parameters  of  a  product  are  analytic  functions  of  the 
parameters  of  each  factor  in  the  product.  For  example. 


L(y)  =  L(a)L(P) 


where 


y  =  f(a,p). 

This  guarantees  that  the  group  is  differentiable.  The  Lie  groups  used  in  Yang-Mills  theory 
are  compact  groups,  i.e.,  the  parameters  range  over  a  closed  interval. 


13.1  Calculations: 

The  nonlinear  Schrodinger  equation  (NLS)  is: 


where 


q  is  the  normalized  effective  amplitude  of  the  wave  electric  field: 

*-2a**. 


\ 

J 


where  g  is  the  gain,  is  the  Kerr  coefficient,  A  =  2 ncla,  e-  Lco/cOy  is  a  small  parameter 
and  E  is  the  electric  field  intensity; 
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T  is  time  normalized  by  the  characteristic  time,  of  the  side  band  spectrum,  (AtOj)  \and 

the  coordinate  %  defined:  f  =  e(t  —  k'  z),  where  k’  is  the  derivative  of  the  wavenumber  or 

dk/dco  evaluated  at  cop  the  angular  frequency  of  the  coherent  light  carrier  and  z  is  the  axial 
distance.  T  is  therefore: 


T  =  ~. 


u 

Zis  the  axial  distance  normalized  by  the  dispersion  distance,  z0  =  -t2Jk",  defined  for  the 
spectrum,  A cox,  and,  £,  defined:  g  =  e2jz.  Z is  therefore: 


Z  = 


_  g 


A  =  is  ^  linear  hgher-order  dispersion  coefficient, 

o  tQ 

The  bound  two-soliton  solution  or  A/- soli  ton  solution  of  the  NLS  equation  is3: 


q  = 


where 


2fZL±ZkA 

-HtJ 


{’ll  cosh  e2  exp[zo\]  +  t)2  cosh  6X  exp[za2]) 


A  =  cosh(0j  +Q2)  + 


rIk± Jh} 

Vl  ~~V2  J 


cosh(9,  -  e2)+ -iaa_cos(cri  _  2), 


dt  =  i i,T+el0  i  =  1,2 
90l  =  constant. 


The  phase,  a„  satisfies: 


da,  1  , 


1 1  represents  the  amplitude  as  well  as  the  pulse  width  of  the  soliton. 

This  solution  is  found  for  the  input  pulse  shape  initial  value: 

q(T,0)  =  A  sec  hT. 

Fig.  13.1  shows  examples  for  two  values  of  rj1  and  i i2  and  with  9l  =  ^T,  at  =  rfz /2  for  / 

=  1,2. 


J  ^  Higher-order  terms  appear  at  distances  longer  than  1  &z0.  In  this  case,  77, ,  are 
defined  as3:  u 


3  Hasegawa,  A.  &  Kodama,  Y.,  Solitons  in  Optical  Communications,  Oxford  U.  Press,  1995. 
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7?i>2  =  l±2exp  ^ 

Sspersionh°WS  ^  mteraction  of  two  ^ual  a“Plitude  Pulses  in  the  presence  of  third-order 

t  ^  joider  to  ^crease  the  rate  of  pulse  generation  it  is  desirable  to  place  solitons  close 

°  eacl;  ®thfr-  However,  because  of  the  nonlinearity  of  the  system,  solitons  mutually 
interact,  leading  to  deterioration  in  the  repetition  rate.  In  the  case  of  in-phase  solitons  with 
equal  amplitude  we  consider  the  solution  of  the  NLS  equation  with  two  soliton  interaction! 

w«h  T0  =  7,  corresponding  to  ry1  =  1.072,  T)2  =  0.952  and  ZD  =  26.3  (Fig.  13  3)  Z  =  Z  /2 
is  the  collision  distance  of  the  pulses.  ^  8  '  c  Zp/Z 

asymptotic  °f  °ff'phase  Solitons  ™th  ^  amplitudes,  the  NLS  takes  an 


2 

q(J,Z)  -  X *7/  sec hT}t(T  +  KtZ  -  J^expf-i^r  + 1(77? 


«?)/2  +  i<4],  cwZ-»±oo. 


where 


v12  =  ±2exp^jsin^J, 
t?ii2  *  1  ±  2exP^- yjcosj^-J 

Fig.  13.4A  shows  the  evolution  of  two  unity  amplitude  seeh  pulses  with  an  initial  phase 
difference  of  4>0  =  nIA  and  an  initial  pulse  separation  of  T0  =  7. 

In  the  case  of  in-phase  solitons  with  unequal  amplitudes. 


*7l,2 


Fig.  i  3.5  shows  the  eigenvalues  rjl  2  for  two  unequal  amplitude  sech  pulses  as  a  function 
of  the  initial  separation  T0  and  Fig.  13.4B.  shows  soliton  interaction  with  two  uneaual 
amplitude  pulses,  A0  =  1.1  and  the  initial  pulse  separation  T0  =  7.  ^ 
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lg.  13.1.Two  examples  of  bound  two-solitons.  A:  t/j  =  3,  rj2=  l,  giving  Zp=  ?t/2;  B:  Tfc  =  1.5,  n,=  1 

giving  Z,  =  1671/5.  12  ’ 


Fig.  13.2.  Interaction  of  two  equal  amplitude  pulses  in  the  presence  of  third-order  dispersion.  (A):  initial 

separation  T0  =  10;  (B)  initial  separation  T0  =  7. 


Soitonintemcioji withttro  equal  turpMuAe pu2ses; 

initial  pulse  nparaiion:  T„  =  7,  correspov&xg  to :  m  =  1.072;  rf2  =  0352;  2,  =  263. 


\q\ 


L. 


Fig.  13.3.  Soliton  interaction  with  two  equal  amplitude  pulses.  Initial  pulse  separation  T0  = 

Evob*ion.»f  Two  Utdfy  AvfUkUU SKk.Pals*:wXkmitiml  jpkmx  Stffvyneeof:  wXkiwo  uaKpml  mmftHmlt  jm2m£ 

<Po  =*/*;  auiml pmit# stpondmx:  T0  -7.  *o  T0  =  7. 


7. 


Fig.  13.4.  (A)  Evolution  of  two  unity  amplitude  sech  pulses  with  an  initial  phase  difference  of  <t>0  =  tt/4. 
Initial  phase  separation  ,  I0  =  7.  (B)  Soliton  interaction  with  two  unequal  amplitude  pulses,  A0  =  1.1  and 
the  initial  pulse  separation  T0  =  7. 
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14.0  The  Bootstrap 


This  section  addresses  analyses  offered  by  the  Bootstrap  method  for  estimating  the 
accuracy  of  a  paramater  estimator  taken,  in  the  instances  addressed,  as  radar  signal  returns. 
Here,  we  show  -  for  the  first  time  to  our  knowledge  -  that  the  technique  can  be  applied  to 
distinguish  targets.  We  also  apply  techniques  taken  from  fractal  theory  and  deterministic 
chaos  theory  to  show  that  radar  returns  have  both  a  fractal  nature  and  are  not  time  series 
produced  by  a  random  mechanism,  but  rather  by  a  deterministic  chaotic  mechanism.  We  are 
able  to  show  that  the  fractal  properties  of  the  power  spectra,  the  variance  and  Fano  Factor 
distinguish  between  two  target  returns.  Use  of  embedding  dimensional  analysis  also 
reveals  the  deterministic  chaotic,  as  opposed  to  random,  nature  of  the  radar  returns.  The 
third  report  from  Northwestern  University  is  also  included. 

The  bootstrap4  is  a  technique  for  estimating  the  accuracy  of  a  parameter  estimator 
where  additional  re-samplings  to  gain  additional  raw  data  is  not  possible.  In  radar  signal 
processing,  large  sample  methods  are  not  possible.  Here  we  apply  this  technique  to 
determine  differences  between  two  target  returns  (Fig.  14.1).  The  technique  has  only 
commenced  to  be  used  in  radar  processing5.  The  bootstrap  does  what  in  an  ideal  situation  a 
radar  operator  would  do  if  it  were  possible  -  that  is,  the  operator  would  repeat  again  and 
again  irradiating  a  stationary  target  (as  if  it  were  on  a  range).  In  the  case  of  the  bootstrap, 
the  measurements  in  each  time  bin  of  the  signal  return  are  randomly  reassigned  and 
estimates  are  recomputed.  This  can  be  done  hundreds,  even  thousands,  of  times  and  the 
results  are  treated  as  if  they  were  real  repeat  test  measurements. 

In  Ref5  there  is  discussed  an  application  in  which  the  distribution  of  the  estimated 
“close  approach  probability”  is  used  as  an  index  of  collision  risk  in  air  traffic  control.  Here, 
we  show  -  for  the  first  time  to  our  knowledge  -  that  the  technique  can  be  applied  to 
distinguish  targets. 


4  Tibshirani,  R  J.,  Variance  stabilization  and  the  bootstrap.  Biometrika,  75, 433-444,  1988. 

Efron,  B.  &  Tibshirani,  R.,  Bootstrap  methods  for  standard  errors,  confidence  intervals  and  other  measures 
of  statistical  accuracy.  Statistical  Science,  1,  No  1, 1217-1241, 1989. 

Efron,  B.  &  Tibshirani,  R.J.,  An  Introduction  to  the  Bootstrap,  Chapman  &  Hall,  NY,  1993. 

Zouhir,  A.M.  &  Boashash,  B.,  The  bootstrap  and  its  application  in  signal  processing.  IEEE  Signal 
Processing  Magazine,  January,  15,  No.  1, 56-76, 1998; 

Politis,  D.N.,  Computer-intensive  methods  in  statistical  analysis.  IEEE  Signal  Processing  Magazine 
January,  15,  No.  1, 39-55,  1998;  "  " 

5  Nagaoka,  S.  &  Amai,  O.,  A  method  for  establishing  a  separation  in  air  traffic  control  using  a  radar 
estimation  accuracy  of  close  approach  probability.  /.  Japan  Ins.  Navigation  82, 53-60, 1990; 

Nagaoka,  S.  &  Amai,  O.,  Estimation  accuracy  of  close  approach  probability  for  establishing  a  radar 
separation  minimum.  J.  Navigation  44, 110-121, 1991. 
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Fig.  14.1.  Two  radar  return  signals,  A.  and  B. 


The  basic  principle  is  that  if  X  =  {X  j,X2,...,Xn}  is  a  sample  -  in  the  present 
instance  it  is  a  radar  return  signal  -  it  may  be  considered  as  a  collection  of  n  numbers  drawn 

at  random  from  an  unspecified  distribution,  F.  If  denotes  a  characteristic  of  F,  e.g.,  its 

mean  or  variance,  then  if  an  estimator  of  were  known,  then  ri  could  be  judged  to 

exceed  a  certain  bound. 

Put  another  way,  the  boostrap  procedure  assumes  that  the  sample  X  itself  is  die 
underlying  distribution.  By  resampling  from  X  many  times  and  computing  t?*  for  each  of 
the  resamples  a  bootstrap  dsitribution  of  is  obtained  that  approximates  the  real 
distribution  of  ri*  and  from  which  confidence  intervals  can  be  derived. 

A  way  this  can  be  done  is  if  the  X’s  are  considered  to  be  independent  and 
identically  distributed  random  variables  and  are  part  of  the  sample  distribution,  F*,  which 

is  considered  close  to  the  true  distribution,  F,  i.e.,  F*  approaches  F  as 

Using  a  pseudorandom  number  generator,  100  samples  of  177  values  (Signal  #1) 
and  812  values  (Signal  #2)  were  drawn,  with  replacement,  from  the  Signals  #1  and  #2 
(Fig.  14.1).  The  means  and  standard  deviations  of  these  samples  were  then  calculated  and 
sorted.  The  results  are  shown  in  Fig.  14.2,  with  mean  standard  deviations  and  confidence 
limits  indicated. 
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Fig.  14.2.  Bootstrap  Standard  Deviations:  100  random  samples,  with  replacement,  from  the  radar  signal 
returns.  In  A.  the  mean  standard  deviation  is  0.0476  and  the  95%  confidence  interval  is  (0.0419, 0.0539).  In 
B.  the  mean  standard  deviation  is  0.1154  and  the  95%  confidence  interval  is  (0.09553, 0.15114). 
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The  bootstrap  method  also  is  able  to  indicate  skewness  of  distributions.  Fig.  14.3 
shows  histogram  plots  of  the  data  shown  in  Fig.  14.2  with  superimposed  probability 
density  function  of  a  Gaussian  variable.  It  is  evident  that  Signal  #1  is  skewed  to  the  higher 
density  values  and  Signal  #2  to  the  lower. 
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Fig.  14.3.  Histograms  of  the  bootstrap  standard  deviations  based  on  random  samplings  of  signal  A.  and 
signal  B.  The  solid  line  indicates  the  probability  density  function  of  a  Gaussian  variable  with  a  mean  of  5.0 
and  variance  2.5. 
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14.1  Fractal  Dimensions  and  Chaos 

The  power  spectra  of  fractal  processes  in  time  reflect  the  statistical  properties  of 
those  properties.  Self-similarity  means  that  there  is  a  relationship  between  the  power  at  the 
high  frequencies  (fine  resolution)  and  the  power  at  the  low  frequencies  (coarse  resolution). 
This  relationship  has  power  law  scaling  and  the  energy  at  a  given  frequency  is  proportional 

to  1  If.  Therefore,  log( power  spectrum)  =  -a  fog(spatial  frequency),  and  a  can  be  found 
from  the  slope  of  the  power  spectrum  curve  in  log-log  plots. 

In  the  case  of  fractal  processes,  the  total  power  at  all  frequencies  in  the  power 

spectrum  and  the  average  power  do  not  exist.  If  a  >  1,  (as  in  Fig.  14.4,  below),  the  total 
energy  increases  as  the  lowest  frequency  used  to  measure  it  decreases,  or,  the  longer  the 

interval  of  data  that  is  analyzed,  the  larger  the  total  power  in  the  power  spectrum.  If  a  <  1 , 
the  total  power  increases  as  the  highest  frequency  used  to  measure  it  increases,  or,  the  ever 
shorter  intervals  of  data  contain  ever  larger  amounts  of  power. 
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Fig.  14.4.  Log-Log  plots  of  spatial  frequency  versus  the  power  spectrum.  The  solid  line  indicates  the  fractal 
power  spectrum  relation:  P(f)  =  Ilf.  The  slope  of  A  is  -2.98985;  of  B  is  -2.89493. 

In  the  case  of  fractal  processes,  the  moments  (mean,  variance,  etc.)  depend  on  the 
U?ed  t0  measure  th®m-  Therefore,  multiple  sampling  of  the  signals  shown  in  Fig 

in  twf  enraVei^g!.temp0r^  samplmg  ^dows,  reveals  the  dependent  of  the  moments 
on  the  sampling  window  resolution.  The  scaling  relationship  describes  this  dependencv  on 

Jnr^nf  fr  1S  cKharactenze/  by  *e  fractal  dimension,  d.  Again,  ctticulatuJ  die 
slope  of  a  log-log  plot,  this  time  of  the  moments,  which  for  Fig.  14  5  are  the  standard 

versus  *e  sampling  window,  provides  the  method.  The  SeriS 

values  of  d  are  given  in  the  caption  to  Fig.  14.5.  8 


~  - -  vcisus  sampling  nme  windows.  The  solid  line  indicates 

how  a  moment  -  in  this  case  the  standard  deviation,  depends  on  the  resolution  (temporal  sampling  window 
duration)  used  to  measure  it.  The  slope  of  A  is  0. 1385;  of  B  is  0.05425.  S 

According  to  Takens’  theorem6,  the  entire  phase  space  set  of  a  deterministic  chaotic 
contracted  from  one  independent  variable.  A  consequence  of  the  theorem  is 
that  a  deterministic  chaotic  process  will  exhibit  a  phase  space  set  with  increasing  N  the 
dimension  of  the  phase  space,  or  the  number  of  values  of  the  data  taken  at  a  time  tefproduce 
die  phase  space  set  (the  embedding  dimension ).  Thus,  the  S  dtosTon  oH 
deterministic  chaotic  process  will  reach  a  limiting  value.  However,  if  the  process  is 
random,  the  fractal  dimension  of  the  phase  space  set  increases  as  the  embedding  dimension 
increases,  i.e.,  as  the  number  of  values  of  the  data  taken  at  a  time  to  produce  the  phase 


3^198’/"  Rand  &  Y°Ung>  (Ed‘S)  °ynamical  Systems  Turbulence,  Springer-Verlag,  NY,  pp.  366- 
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SnSaS  JfCTeases‘  1x1  cfse  of  **  deterministic  chaos  process,  the  time  series  was 
by  a  fmite  ^  °f  ^dependent  variables,  and  in  the  case  of  a  random  process  the 
time  senes  is  generated  by  an  infinite  set  of  independent  variables.  process,  the 

c.  w  j4-6  shows  that  as  the  embedding  dimension  increases,  both  signals  do  not 

slishtiv  lestoi  dimentsi0"'  ?a,her-  bodi  signals  reach  Ifaitfag  tin™ 

mif  Z  ^  k  ^  SJ  Ume  senes  for  both  signals  were  generated  by  a  deterministic 
rule  that  can  be  descnbed  by  3  equations  with  3  independent  variables. 
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Fig.  14.6.  The  fractal  dimension  of  a  phase  space  set  of  a  random  time  series  increases  as  the  emhedrtW 
chmension  increases.  Here,  it  is  shown  that  as  the  embedding  dimension  increases  (=  number  of  values  of 

tUne1t°  S**  1116  Phase  ^  «0.  the  fractal  dimension  if Tphai  Lte  * 
appoachesahmmng  value.  The  fractal  dimensions  are  of  the  power  spectra  as  in  Fig  14  4^towiS 

*•- wi,h  Fis- i4A  ^  or 


,  , As  *e  statistical  properties  (e.g.,  the  mean  and  variance)  of  a  fractal  time  series 

Srisri™?n  thC  reS0  Utl0nn  used  t0  mea?ure  them.  «  is  not  appropriate  to  measure  these 
tistical  properties  at  only  one  resolution.  However,  it  is  appropriate  tn  Hewmin^ 

kef  ProPerties  depend  on  the  resolution  used  to  measure  them.  The  Fano  Factor 

is  a  statistical  measure  appropriate  for  measuring  the  statistical  DroDerries  r»f 
processes  and  is  defined  as  the  variance  divided  bv  the  mean  The  Fann  Farmr  v,0  ^actal 
taw  fonn  that  is  proportional  to  r»,  where"  ? £  SSi  ™d  t^flre  S 

ve?urth?OT?raIta1^?“r?,ffrt0hm  *e  S!0pe  of  ato«to*  plot  of  the  Fano  Factor 
shows  snchSr  SVe  ^ofKlT  ^  °f  5eries'  F*  14  7 

The  fiactal  dimension,  d,  is  related  to  b  by: 


d  = 


_t± 

i 

a 


where  a  and  fi  are  exponents  determining  how  the  measured  property  depends  on  the 

of  and  the  resolution.  We  have  established  (Fig.  146)  tht£  thlprocS 

gener^ng  the  srgnals  is  deterministic  chaotic  and  described  by  3  independenfvSs 

Therefore  fi-3  (and  a  -  1).  From  the  values  for  b  calculated  from  Fig.  14.7  and  shown  in 

ffifS “v  *£££*?  *P  ’  3:09  (Signal  #2),  which  valuTaS 

the  range  of  acceptable^enor!11  t^°Se  va^ues  ^  ^  shown  m  Fig.s  14.4  and  14.6,  or  within 
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Fig.  14.7.  The  Fano  factor,  F,  is  equal  to  the  variance  divided  by  the  mean.  Here  are  Loe-Lno  nintc  «f  a 
Fano  Factor  as  a  function  of  the  temporal  length  of  the  samnling  winHrm/c  tj,  ^  -  ^  ^ 
reladonship  of  the  Fano  factor  with  respect  to  theremporal  sampling  Ldow 4  P°”  A 

*■—  >-)■  »  *  c^of  A.  L  KSS  Wi  •’it"**." 
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15.0  Summary  Recommendations 


As  was  stated  in  1.3  State  of  the  Art  Benchmark  from  which  BSEI 
Research  Proceeds,  we  have  sought  to  either  identify  or  create  methods  different  from 
state-of-the-art  for  solving  the  BMDO  problem  of  providing  an  ultrahigh  resolution  sensor 
which  can  identify  missile  from  decoy.  We  advocate  sensor/radar  systems  using  pulses 
shorter  than  the  length  of  the  target,  adaptive  time-frequency  analysis  method  which 
SSdecoy^§et  backscattenng  centers  and  resonances,  and  pattern  identification  of  target 

In  particular,  we  have  created 


•  extensions  of  the  Gabor  transform  (CQOW) 

•  adaptive  methods 

•  methods  involving  the  instantaneous  phase 

•  mutual  information  methods 

•  Zak  transform  methods 


Furthermore,  we  have  created  methods  using 


•  the  bootstrap 

•  fractal  dimension 

•  Fano  number 

•  deterministic  chaos  methodologies 


to  provide  precision  in  identification  of  target  and  decoy. 


Our  recommendation  for  the  next  step  to  be  taken 
methods  in  emulation  or  real-time  systems. 


is  to  implement  these  new 
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